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qn the other hand, the radiative coupling may 
be increased by as much as a factor of 5 com- 
pared to surface bursts if there 'is a slight pene- 
tration of the weapon into the surface of the 
ground (3 to 5 meters for a I Mt burst), This 
greater coupling is in addition'to the well known 
increase in mechanical coupling, which is dis- 
cussed in paragraph 247. 


2-47 
Mechanical Coupling 
The preceding d'iscussion outlined the flow of 
energy into the ground for cratering events that 
result from bursts at or above the surface. These 
processes take place during the flI'St several milli- 
seconds of actual crater formation. The hydro- 
dyn8JTIic portkn of crater fonnation may last 
for several seconds. 
.. 
. 


'. If the explosion occurs at a height-of-burst 
greater. than 4 to 5 m/kt 1/3, the crater will be 
formed 
almost 
exclusively 
by 
compaction. 
Scouring, which results from airblast an'd airblast 
associated winds, contributes a small fraction of 
the crater volume. Since ther.e is little fallback in 
such a crater, the true and apparent craters are 
essentially identical. As the height of burst is de- 
creased, an increasing fraction of the crater is 
formed by excavation and ejection of ground 
material. 
In 
addition, 
compaction, 
plastic 
deformation and flow of the material in the rup- 
ture zone are principal mechanisms for crater 
furmation 
of very neaNurface bursts. These' 
mechanisms produc'e a slight increase in material 
density as well as considerable brittle and plastic 
failure of the ground material. The me.chanisms 
also cause a radial displacement with 'an over- 
turning of crater material around the edge of the 
crater. The top of this mound is tenned the 
crater lip, After detonation, the crater maybe 
partially remled by material falling back into the 
crater (fallback), by slumping, and by late-time 
rebound. The material thrown o,ut beyond the 
crater lip (ejecta) fonns a layer whose thickness 
diminishes with radial distance from the burst 
point. 


The first diagram in Figure 2-70b illustrates 
the general crater shape that is expected from an 
above-surface 
burst 
Despite many 
explosive 
cratering experiments with. nuclear and high- 
explosive (HE) devices, the complex interrela- 
tionshipsthat exist among the various mecha- 
nisms 
that 
could 
cause .craters preclude the 
quantitative dynamics of crater formation by 
near-surface bursts from being well understood. 


The size and shape of craters from subsurface 
detonations 
are 
affected significantly by the 
depth of burial. The dominant cratering mecha- 
nisms, which are discussed in succeeding para- 
graphs, will be related to the depths of burial 
shown in Figure 2-70b. 


• 
Crushing, Compaction, Plastic Deforma- 
tion. As the high J'ressure explosion gases 
expand against the. surrounding medium, 
a spherical shock wave is generated. The 
initial shock pressUres are as large as 10 to 
100 
million 
atmospheres. The medium 
surrounding a nuclear explosion is ini- 
tially melted and vaporized as the shock 
passes through 
it. 
As the shock front 
moves outward in a spherically diverging· 
shell, the medium behind the shock front 
is put into radial compression and tan- 
gential tension. The peak pressure in the-- 
shock front attenuates by spherical diver- 
gence and dissipative mechanisms in the 
medium. The material around the explo- 
sive charge (nuclear or HE) is crushed, 
heated, and physically displaced outward, 
forming a cavity. In regions outside this 
cavity. the shock wave will produce per- 
manent deformation by plastic now. This 
mechanism of crater fonnation isa sig- 
nificant 
contributor to crater size 
for 
bursts at or just ,below the ground surface 
(FigUre 2·70b,(b»; and at depths of burial 
that result in subsidence craters. (Figure 
2-70b,(0). 
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• 
Spalling. When an upward moving shock' 
(cortpressive) 
front encounters the air- 
ground interface, the large mismatch of ,. 
material properties results in the genera- 
tion 
of 
a 
n.egative 
stress 
(rarefaction) 


wave. 
The 
refraction 
wave 
propagates 
back into the medium in which the burst 
occurred, and puts the medium 
(origin~ 


ally, under high compression) into tension 
(or [less compression). This phenomenon 
causes the medium to break up and move 
up"";ard with a velocity characteristic of 
the total momentum imparted to it. In a 
loose soil material, this spalling ejects in- 
dividual particles of siJ:Tlilar siz" .into the 
air. In a rock medium, however, the thick- 
ness and size of the spalled material gen- 
erally is determined by the presence of 


pre~xisting fracture patterns and zones of 
weakness. The velocity of the spalled ma- 
terial decreases in .proportion to the de- 
crease in 
peak compressive stress. The 
spall 
mechanism produces an extended 
rupture and plastic zone· near the ground 
surface and contributes significantly to 
the true lip height of' the crater. This 


~echanism appears to be dominant in de- 


. termining crater size at shallow depths of 
burial (Figure 2-70b,(c)). 


• 
Gas A cceleration. Gases are produced in 
the material surrounding the explosion by 
vaporization 
and 
chemical 
changes in- 
duced 
by heat and pressure. The near- 
adiabatic expansion of the gases imparts 
motion 
to 
the 
medium. At depths .of 
burial 
at 
which 
crater dimensions are 
maximized (optimum DOB), the cavity 
gases produce appreciable acceleration in· 
overlaying material before 
they escape 
(vent) through cracks extending from the 
cavity to the surface. Cas acceleration is 
the dominant crater producing mechanism 
at optimum DOB (Figure 2-70b,(d)). At 
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shallow depths of burial the gases cannot 
exert significant pressure before ventmg 
occurs. In the case of vel)' deep explo- 
sions, the weight of the overburden pre- 
cludes any significant gas induced accel- 
eration of the overlying material. 


• 
Overburden Collapse. At depths of burial 
that are large (two or more times) com- 
pared to the optimum, the mechanism of 
overburden collapse (subsidence) becomes 
dominant. This effect is closely linked to 
the crushing, compaction and plastic de- 
formatio'n mechanism that produces an 
underground cavity. At these large depths 
of burial, spall and 
~as accelerqtion will 
not impart sufficient velocity to the over- 
lying material to eject it physically from 
the crater. It would' be expected that the 
crater 
volume 
would 
be 
determined 
largely by the underground cavity formed 
by the detonation. In a rock medium, 
however, material that is fractured and 
displaced from its original position tends 
to take up more space than it occupied in 
its natural state. This bulking action could 
result 
in 
no crater or, indeed, even a 


mound above the ground. A mound was 
produced after the SULKY detonation, as 
shown in Figure 2-72a. 
In 
certain 
geologic 
materials 
yet 
another type. of subsidence occurs. When 
the pressure in the cavity decreases below 
overburden 
pressure, 
the 
roof of the 
cavity begins to collapse. In most media, 
this collapse will continue upward and 
will form' a "chimney" of collapsed ma- 
terial. In a soil medium, where the density 
of the material will not change signifi· 
cantly after it has fallen, the volume of 
the cavity will be transmitted to the sur- 
face, forming a "subsidence crater" on the 
surface (Figure 2-70b,(f». This phenom-. 
enon is illustrated in Figure 2-72b. 
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Figure 2-72a. SULKX Event; Mound Created by the Bulking of ROckIMaterial in e Deeply Buried 
0.087 kt Nuclear Detonation: 
Mound Oiemeter" 48.8 m 1160 ft.); Mound Heiaht .. 7.6 m (25 ft.); OOB .. 27 m(90 ft.) 
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Figure 2-72b. PASSAIC Event; Subsidence Crlter Produced by I 
Deeply Buried Nuc\llr Detonetion 
in '8 Soil Medium 
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PREDICTION OF CRATER DIMENSIONS 


The 
present 
prediction 
of crater 
dimen- 
sions for ab6ve-surface, surface, and shallow- 
buried burSts is based'on the nuclear tests con- 
ducted in dry alluvial soil at the Nevada Test 
Site 
and 
the 
saturated 
coral at the Pacific 
Proving Ground. Extrapolation to other geolo- 
gies is made ,on the basis of HE spherical<harge 
events, 
There 
have 
been no high-yield near- 
surface nuclear experiments in dry soil or rock. 
It is assumed that the influence of geology can 
be separated from that of the energy source, 
However, the nuclear tests conducted at the 
Pacific Proving Ground (saturated coral) and HE 
tests conducted at the Suffield Experiment Sta- 
tion (wet soil) suggest that the influence of 
geology 
is 
not 
entirely 
independent 
of the 
source. This fact is illustrated in Figure 2-73. 
The dimensions of the crater vary in· a com- 
plex manner as the HOB approaches the ground 
surface; but 
with 
lower 
HOBs the apparent 
depth of the crater generally increases and the 
ratio of the apparent radius to depth decreases. 
When 
crater 
dimensions 
for 
surface, 
above- 
surface, and shallow-buried bursts are scaled to 
yields much greater than I kt (for example, I 
Mt), the calculated change in the shape of the 
apparent crater as the yield increases is also an 
unportant 
consideration. 
This 
phenomenon, 
unfortunately, has not been observed in dry 
media, 


248, Crater Dimension Scaling 


The results of cratering explosions of various 
yields and burst geometries have been corre- 
lated 
by 
empirical' scaling' laws that express 
crater dimensions in terms of a standard yield of 
one kiloton. Crater dimensions for any other 
yield can be predicted by application of an em~ 
pirical, yield scaling exponent, as will be de- 
scribed in the succeeding paragraphs. The dimen- 


sions for which curves and scaling procedures 
will be presented ,are the radius of the apparent 
crater, R a , the depth of the apparent crater, 


,Da, and the apparent ·volume ~. Other pertinent 
crater dimensions such as the radius to the crest 
of the apparent lip and the height of the appar- 
ent lip, may be related to the above dimensions 
by scalmg laws that will be described. 


Several burst geometries will be considered: 
above-surface (near and contact), surface, shallow- 
buried, 
and 
deep-buried 
bursts. 
The 
actuai 
height of burst (HOB) or depth of burst (DOB) 
is measured from the original ground surface to 
\ 
' 
the center of energy of the weapon. The scaled 
HOB' or DOB is the actual dimension divided by 
the scaled yield (described belo'.,,). Of the above- 
surface bursts (HOB/Wl!3 > 0), the lower bursts 
are 
of 
most 
significance 
to 
cratering, 
i.e., 
HOB/WI!) 
,<;; 3 m/kt1/3. A contact burst is 
considered to be one wherein the weapon is in 
contact with the ground surface, and therefore 
has an actual HOB of approximately 0.5 m. A 
surface burst is one in which the HOB as defined 
above, is exactly zero.' A shallow-buried burst 
has a DOB/Wl!3 <:; 5m/kt1!3 and the deep-buried 
burst has a DOB/WI/) > 5 m/ktl/3. 


No single yield scaling exponent has been 
found 
to be valid for scaling apparent crater 
dimensions over a wide range of yield, geology, 
and 
HOB or' DOB. 
However, 
the 
following 
yield scaling exponents, a., (yield, W(kt), to the 
a. power) have been found to be approximately 
correct for 
scalin~ the apparent crater dimen- 
sions for' surface/above-surface bursts and for 
deep-buried bu~ts: 


ci= 1/3 for 0 <; HOB <; 3 rn/kt 1!3 


kllI3 
.' 


a. = 1/3.4 for~ > 5 m/ktl/J• 
.....1!3 
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Figure 2-73. 
Typical Pacific Proving Ground Crater Profiles 
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in the shallow-buried region where the DOB/WI/3 


is less than 5 m/kt l!3, the yield scaling exponent, 
n, is a function of depth of burial.· It is in this 
region that 
1an interpolation procedure is re- 
quired to determine a, as will be discussed in 
paragraph 2-49. 
The 
near-surface 
region 
(HOB/WI/3 
EO;; 
3 
m/kt l/3 to DOB/WI/3 
EO;; 5 m/kt 1/3 ) and th'e 
buried burst (DOB/Wl/3 > 5 m/kt l /3) region 
will be examined separately in the succeeding 
paragraphs. 
Five 'generic, geologies have been 
chosen to describe the effect of various homo- 
geneous media on cratering efficiency. 


2..4.9 
Volur-e of Craters from Surface 
and Near-Surface Bursts 
(HOBIW 1/3 ~ 3 m/kt1/3 to 
DOB/W1/3 ~ 5 m/kt1f3 ) 


Figures 
2-74a through 
2-75e show the ap- 
parent crater volume, as a function of height 
of burst,t for a near-surface I kt explosion in 
various homogeneous generic geologies. Figures 
2-74a through 
2-74e are valid for yields less 
than or equal to I kt and Figures 2-7 Sa through, 
2-75e are valid for yields greater than 10 kt. 


For yield s between 
I and 10 kt, the crater 
volume may 
be estimated 
by using the best 


p-~tim;jte curve of the appropriate geology for 
either I kt or 10 kt or a weighted average, de- 
pending upon the output characteristics of the 
warhead. The higher the radiative fraction, the' 


·Craters from detonations on the islands of the Pacific Proving 
Ground (Eni\l'etuk and Bikini atolls) indicate that, for a aatu- 
rated 
coral medium, the scaling exponent for the apparent 
criter radius may be greater than 1/3 and for the crater depth, 
smlller than ,1/3. Until the reason that it is necessary to depan 
from cube-root scaling in aaturated coral is understood. it will 
not be possible to determine £imilar depanures in other media, 


t in this group of fl8ures. HOB is used to indicate both HOB and 


~.:;r., ..e., buried bursts are shown with a negative HOB. 


" 


closer the crater volume will be to the W> 10 kt 
curve. If, for example, it is known that a 3 kt 
weapon of interest has a high radiative output, 
the W > 10 kt curve should be used to deter- 
mine the apparent crater volume. If, however, 
no information is available concerning the radia- 
tive output, the crater volume should be deter- 
mined by interpolation, assigning a 60 percent 
weight to the W <; 1 kt curve and 40 percent 
weight to the W > 10 kt curve, i.e., for a 3 kt 
weapon with no further information 


y'(3) =0.6Va(W< 1)+ 0.4 Va(W> 10). 


The 
yield 
scaling parameter, a, 
for 
near- 
surface bursts is determined as follows: 


HOB;> 
0·0. = 0.33'3 
Wl/3' 
' 


HOB 
-5 ~- < 0 (m/kt 1/3 ). a =(see Figure "'I-76a) 
k!1/3 
' 
,.' 


Once the value of a is known, the apparent 
scaled 
crater 
volume 
can 
be 
obtained 
from 
Figures 2-74a through 2-75e for a specified yield 
(W) and HOB/Wa 
in the appropriate geology. 
The value, obtained by 
this procedure is the 
"best estimate" of the apparent crater volume 
for a I kt explosion (¥al)' 
' 
The crater volume for yields other than I kt 
may 
be 
obtained 
by' 
the 
following 
scaling 


relation; 
v' 
_a = W3a 
Y.I 
• 


where Val is the apparent crater volume for a I 
kt explosion, and V. is the corresponding vol- 
, ume'for a yield of W kt. The height of burst 
scaling relation is 


HOB 
HOB =wa. 
1 


where HOB 1 is the height of burst for a I kt 
explosion, and HOB is the corresponding height 
of burst for a yield of W kt. 
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The u:ncertainty in the crater volume for W <; 


I 
Ie.t is I based on HE criltering results. It is 
assumed that 
theun~rtaintyis independent of 
geology :and yield. The bounds shown in Figures 
2-74" through 2-74e are not the results of a de- 
tailed st,atistical analysis. These bounds contain 
95 percent of the experimental data, but statis- 
tical implications should not be drawn from this 
fact. 
The best estimate for the W> 10 kt curve was 


det~rmined by nuclear data from the Pacific .' 
Proving 'Ground. The lower bound was suggested 
by calculations. and was influenced by"the fact 
that all the nuclear data were obtained from de- 
,ices 
t]:lat 
were 
very 
dissimilar 
to 
modem 
weapons. 


2·50 
Crater Shapes and Dimensions for 
S!Jrface and near·Surface Bursts 


The Shape of the crater from a near-surface 
burst must be known to determine the radius 
and depth from the crater volume and the ex- 
plosive :yield. The shape depends mainly on the 
yield, the 
~caIed height of burst, and the geol- 
ogy 
Table 2-8 shows the shape to be expected 
as a function of these three variables. The suc- 
ceeding' discussion provides the necessary scal- 
ing relationships to determine the dimensions 
for the various shapes: Separate procedures are 
provIded tor each of the yield/HOB combina- 
tions shown in Table 2-8. 


Table 2-8. Classification of Crater Shapes from 
Near Surface Bursu 85 a Function of 'field, 
Height of Burst, end Geology 


a. W<1 kt 


The crater shape for near-surface bursts (-5 
m/ktl/3 < HOB/Wl/3 <; 3 ,m/ktl/3 ) in all geol- 
ogies for W <;; I kt is best described as bowl- 
shaped. The dimensions of the apparent crater 
radius, R., and apparent crater depth (D.) for 
these 
bursts are 
given 
by the following ex- 
pressions: 


J.J V.1I3 <;; R. <;; 1.4 v.1I3, 


0.35 V.1/3 <. Da ..,;; 0.7 ¥a1/3., 


The best estimates for Ra and Da are 


R a = 1.2 V.1/3, 


Da =0.5 ¥a1/3, 


where V. =the best estimate apparent volume. 
a 
' 


b. W> 1 kt 
The crater shape for near-surface bursts in 
which W> I kt depends on the yield and HOB. 


(]) HOB/k'1/3 ~ O. 
For above-surface bursts (HOB/W1/3 ~ 0), in 
which W > I kt, the crater shape is best de- 
scribed as dish-shaped and the crater dimensions 
are determined by the expressions 
'R 
' 
1.1 WO.OS <; V.U3 <; 1.4 Wo.os, 
a 
D' 


0.35 W-4J·12 ~ v.li3 ";;:0.7 W-OJ2 


• 


All 
All 


HOB/ M 13 ;> 0 
All 


-5 <;; HOB/WI/3 <0 Unsaturated 
Saturated 


Yield 


W<;; I kt 
W> I kt 
W> 1 Itt 
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Scaled HOB 


(m/kt 1/3) 
Geology 
Crater' 
Shape 


, Bowl 


Dish 


Dish/Bowl 
Dish 


The best estimates for Ra and Da are 


R. 
_ 
•. IIl'.OS 


V.1I3 - 1.2",w 
, 
• 
D. 
2 
v.m- =0.5 W-O.l 
. 
• 
where Va =the best estimate apparent volume. 


) 


".,. 
I: 


(2) HOB/W1l3 < 0 


For near-surface buried events 
'(0 < D~! <;5m/ktay, 


in which W >kt, the shape of the crater may be 
either "bowl" or "dish" in nature or some com- 
bination of the two extremes. Mechanisms such 
as compaction, rebound, bulking, slumping and 
iayering effects playa significant role in the fmal 
crater shape. Cratering experience to date is not 
sufficient 
to 
make 
quantitative 
judgments 
regarding the impact of each mechanism. 
The 
apparent 
crater dimensions 
for 
near- 
surface buried 
events are determined by the 
expressions 
Ra 
'Y.'"m =1.2 Wb meters, 
a 
D 
a 


. 


V. 1/3 = 0.5. w-c meters. 
a 


The values of the exponents (b, c), as a function 
of DOB, are found in: Figui'es 2-76b and 2-76c. 
These 
exponents 
are 
necessary 
to 
ensure 
a 
smooth continuous curve for the radii of large 
yield (I 
Mt) bursts as the depth of burial in- 
creases. Existing data indicate that above-surface 
bursts of high yield nuclear explosions produce 
dish-shaped 
craters and buried nuclear explo- 
sions produce 
bowl-shaped 
craters. However, 
since no data exist in the region 


DOB 
0< wa 
<; 5 (m/kto-), 


it is necessary to provide a transition 'region in 
which the crater produced is neither dish-shaped 
nor bowl-shaped, but some intermediate of the 
two 
extremes. 
The above equations prod.uce 


this 
smooth 
transition 
region for the crater 
dimensions. 
It has been noted, however, that in highly 
saturated media slope failure and/or liquefac- 
tion can produce a shallow crater for buried, 


t:~~nts. High explosive tests conducted' at Fort 


Polk, 
Louisiana 
are good examples of these 
phenomena. 
Other 
pertinent crater dimensions may be 
related to the above dimensions as follows: 


• 
The radius to the crater lip crest is 


R aR =J.2SRa· 


• 
The height of the apparent lip is 


.25D. <:'H.~ <:. .33D. 


(near-surface bursts). 


A summary of the procedures for calculating 
craters from near-surface bursts is given on page 
2·161, 
immediately 
preceding 
the 
applicable 
figures. 


2-51 
Dimensions of CraterS from 
Deep·Buried Explosions 


Cratering experiments with chemical explo- 
sives (HE) and nuclear explosions (NE) suggest 
that linear dimensions of craters from buried 
explosions scale. according to a modified over- 
burden rule. This scaling rule is somewhat cum- 
bersome to use for the computation of crater 
dimensions, and the available data do not con- 
clusively prove the validitY of any single scaling 
rule. Therefore, a simplified yield scaling ex- 
ponent, yield (kt) to the 1/3.4 power, has been 
chosen for scaling apparent crater dimensions 
for scaled DaBs greater than 5 m/kt 1/3. Figures 
2-77 through 2:-81 may be used to obtain ap- 
.parent crater volumes for buried bursts in each 
of the soil types indicated. The curves in these 
figures are valid for all yields. The uncertainty 
in crater volume is independent of geology and 
yield, and is based on HE cratering results. 


Crater radius and depth are given by the fol- 
lowing expressions 
. 


1.1 
~113 < R. < 1.4 v.l/3 , 


0.35 ~1/3 <D. <; 0.7 ~113. 
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These 
expressions 
are 
valid 
for deep-buried 


bursts in: all geologies for all yields. The best 
estimates for the crater dimensions are 


R. = 1.2 v.1f3, 


D. =0.5 V.113, 


where v.' is the best estimate of the apparent 
volume, 
It is possible that slumping of the walls will 
produce a wide, shallow crater for some yield- 
DOB combinations in a saturated geology. The 
resulting. crater. shape is neither "bowl" nor 
"dish" in nature; it is some combination of the 
two. The only explosions to date that resulted in 
such occurrences in the craters have been in a 
very wet: clay soil with the explosive charge at 
or near trye optimum depth of burial. 
The 
uncertainty 
in 
crater 
c\;mensions ob- 
tained from the expressions given above is, how- 
ever, approximately correct for any specified 
crater volume. If the crater volume is chosen as 
the 
lowe: 
or upper uncertainty 
band 
from 
Figures 2-77 through 2-81, then the expressions 
for the radius and depth are still valid. A new 
range of crater dimensions is determined by 
using th~ equations and the chosen volume. 
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Figure 2-82 contains the best .estimate predic- 
tion for high-explosive detonations m the five 
gene:ic 
homogeneous geologies. These curves 
show the relative crater volumes for the geolo- 
gies. The curves can be 'used as a guide as to the 
proper nuclear curve to use, for a new or un- 
known geology. If some HE eratering data are 
available, the eurves of Figure 2-82 can be used 
to 
determine which homogeneous geology is 
best for the problem. This figure should not be 
used 
to 
calculate crater volume 
for nuclear 
explosions. 


Other pertinent crater dimensions may be re- 
lated to the above dimensions as follows: 


• 
The radius to the crater lip crest is 


R.~ =1.25R•. 


• 
The height of the apparent lip is 


.20Da <;HIR <; .25Da 
(deep-buried bursts). 


A'summary of the procedures for calculating 
craters from deep-buried bursts is given on page 
2·179, 
immediately 
preceding the applicable 
figures. 


SUMMARY OF PROCEDURES FOR CALCULATING CRATERS FROM NEAR·SURFACE BURSTS 
(-5 m/kt'l'J E;;; HOBIW'f3 E;;; 3 m/kt'!3) 


, 


~' I 


(" 


Figures 2-74a through 2~75e and 2·76a, b and 
c, together with appropriate scaling laws, may be 
used to obtain crater volumes and dimensions 
froll' near-surface bursts in a variety of situa- 
tions. The following summary provides a 
step~ 


by-~t(;p procedure for obtaining such informa- 
tion.. Example problems are provided on the 


page~ immediately follOWing the figures. 


Required (either given or estil?U1ted) 
Yield (H') in kt 
Actual HOB or DOB (may be zero if required) 
Soil type (one of· the five generic types ora 
combination; see Figure 2-86b). 


!. Determine yield 
scaling parameter, a, for 
burst position: 


a. HOB/WI13 ~ 0"" a = 1/3 ~ 0.333, 
b. 0 < DOB/WI/3 E;;; 5 m/kt l/3 .... a from 
Figure 2-76a. 


Actual HOB 
2. Compute HOB I (HOB for I kt) = 
li'Cl 


or DOB~ (DOB for I kt) = ActU~;POB 


3. Determine 
Val 
(apparent-crater volume for 
1 kt): 
a. If 11' E;;; I kt, lise curve for given soil type 
from Figures 2-74a through 2-74e. 
b. If k' > 10 kt, use curve for given soil type 
from Figures 2-75a through 2-75e. 
c. If I kt < Ii' 
E;;; 10 kt, with low or un- 


known radiative output, interpolate using 
Val =0.6 Val (from W <;1 kt curve) + 0.4 


J~ (from W> 10 kt curve). 
d. If I kt < W ,..;; 10 kt with known high 
radiative output, use curve for given soil 
type from Figures 2-75a to 2-75e. 


4. Compute ~ (actual crater volume for yield 
other than I kt): 


V. (actual volume) =Val (W3Cl ). 


These curves are provided for each soil type; 
best 
estimate, 
upper 
bound, 
and 
lower 
bound. Using the following guidelines, select 
the curve most appropriate to the problem. 
a. Use best estil?U1te Val if no specific geo- 


logic data are given. If actual HE tests are 
available for a specific site, compare them 
with the HE cl.:;ves in E.;ure 2-82, and 
move toward the upper or lower bound 
NE curve to adjust for these data, depend- 
ing on where the data are grouped. 
b. Use lower bound 
Val 
(if desirable) for 


target-oriented 
calculations 
(offensive- 
conservative) to calculate R. or Da. Con· 
sider iower-bound Val for targeting when 
calculating crater volume-related phenom- 
ena, such as ejecta, transient velocity or 
displacement. • 


c. Use upper-bound 
Val 
(if desirable) for 


design-oriented calculations (defense-con- 
servative) to calculate R. or Da. Consider 
upper bound Val for design when calcu- 
lating crater-related phenomena affected 
by volume. 


d. For W > 10 kt, the discrepancy between 
the theoretical calculations and the empir- 
ical 
data 
in our understanding of the 
energy 
coupling 
produced 
by 
modem 
weapons detonated at low HOBs or in 
contact· with the earth (HOB = 0.5 m) 


-In all inslaDces it is I'eCOIIIJIleDCl that upper IDd lower bound 
values be calculated in order to pin an appreciation or the 
effects or UDClelUiDties OD the particuJar problem uDder con· 
aideralion. 
. 
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. : 


·See Fi(IIres 2-76b and 2-76c for determining band c. 


must be considered. (See paragraph 2-49.) 
Since this discrepancy only occurs in the 
lower-bound V. values, it principally in- 
fluences targeting, or offensive-eonserva- 
tive', problems. 
. 
It is recommended that the lower-bound 


V. values be considered in calculating kill 
probabilities (Pk ) due to the crater, but 
that 
hi~h credence in cratering Pk not be 
used in assessing overall system Pk . 
Note.' This cratering discrepancy is a sys- 
tematic uncertainty .and should not be 
treated as a random uncertainty. 


5. Calculate R a and Da using the expression 
from the following table that fits the 
}'i~ld 


and HOB and inserting the appropriate v.: 


- 
. 
I 
'- . 


for near-surface bursts.t 


t An exception to tlili is a limilar crater formed by ian explosion 


Oil the Pacific Proving Ground. The crater had no lip and greatly 
reduced ejecta: 11 is thought that such a crater Shape may have 
been the result of a late-limerecoll.lolidation{liquefaction pro- 
cess, which caused the crater lip to be below the waler liurface. 


• 
The radius to the crest of the apparent- 
crater lip is 


• 
The height of the apparent lip is 


The dimensions R a , Da, and Va can be related 
to other pertinent :crater dimensions as follows: 


Variations in the thickness of the ejecta as a 
function of range from the surface ground zero 
(SGZ) are discussed in paragraphs 2-52 and 2-53. 


same 


Below Surface 


(~W~sm/ktI/3) 


Above Surface 
HOB 
(WI/) ;> 0) 


Near-Surface Bursts 


R =I'" Y 1/3 
a 
.6.0 
a 


Da = 0.5 Val/3 
same 


R 
= I 2 WO.08 v. 1/3 
R 
= J.2 wt' V 113 
a' 
a 
a 
a 


D =0 5W~·12V 1/3 D =OSW-c V I/3 
a' 
a 
a' 
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Yield 
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Figure 2-748. Apparent Crater Volume for a '-Kiloton Near-Surface Explosion in Dry Soil; 
Applicable for W __ , kt 
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Figure 2-74b. 
Apparent Crlter Volume for a '·Kiloton Near-5urfece Explosion in 
Wet Soil; Applicable for W";;; , kt 
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Figure 2-74c. 
Apparent Crater Volume for a 1 Kiloton Near·Surfece ExplOlion in 
Dry Soft Rock; AppliCibie tor W< 1 kt 
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Figure 2-74d. 
ApPlrent Crater Vohome for I '·Kiloton Nelr·Surfece Explosion in 
Wet Soft Rock; Appl iCible for W <; , kt 
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( 
, Figur~ 2-74e. 
Apparent Crater Volume for I , Kiloton Nur·SurflCe Explosion in 
Hlrd Rock; Ai!pliCible for W iIi;, kt 
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Figure 2-75a. 
Apparent Crater Volume for a 1 Kiloton,Near-SurflCe Explosion in 
Dry Soil; AppliCible for W > 10 kt 
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Figure 2-7Sb. 
Apparent Crlter Volume for I 1 Kiloton Nllr·Surface Explosion in 
Wet Soil; Applicable for W> 10 kt 
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Figure 2-7Sc. 
ApPirent Crlter Volume forI' Kiloton Nelr·Surfece Explosion in 
Dry Soft Rock; Appliceble for W> 10 kt 
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Figure 2-75d. 
Apparent Cr.ter Volume for .,·Kiloton N••r·SurflCe Explosion in 
Wet Soft Rock; AppliCible for W> 10 kt 
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F.igure 2-75e. 
Apparent Crlter Volume forI' Kiloton Nelr·Surfeee Explosion in 
Hard Rock; Applicable for W> 10 kt 
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Figure 2-761. 
Volume Scaling EJCporent, a. IS. Function of Sclled DOB for 
Near-Surfeee ExplosioN, All Vields 
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Figure 2·76b. 
Yield Transition Exponent. b. for Calculation of the 
Radius of Near·Surface, Buried Explosions W> 1 kt 
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Fi""re 2-76c. 
Yield Transition Exponent, C, for C.lculation of the Depth of 
Nelr-5urfece. Burieo Explosions, W> 1 kt 
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Problem 2·22. 
Crater Dimensions for Above-Surface Bursts 
(HOB >0) 
. 


1'°' 


·1 


.J 


, 
Figures 2·74a through 2-75e include curves 
that indicate the apparent crater volume as a 
function Of HOB for a 1 kt explosion in various 
homogeneous media. 
Scaling. For yields other than I kt, the crater 
volume scales as follows: 


Val =W3a , 


where 
"~' is the apparent crater volume for a I 
kt explosi'on, and Va is the corresponding volume 
for a yield of W kt. The height ofburst scales as 
HOB 
~=wa 
I 
Example: 
Given: ,A hypothetical 300 kt contact burst 
on dry soil. 
Find. The apparent crater radius and depth. 
Solution: The HOB for a contact nudear ex- 
plosion is '0.5 meters. From paragraph 2-49, 
a =0.3333 


for HOB ~ O. Therefore, 


HOB 
- HOB - -.2.:.L - 0 075 
/k If3 
I. - 
Wa 
- (300)lf3 -. 
m 
t 
. 


From Figure 2-75a, the apparent crater volume 
for a I kt explosion at a scaled HOB of 0.075 
m/ktl/3 in dry soil is 


Val =1.25 x 103 m3/kt 3a 


Ya = Val (k')3a c 1. 25 x 103 (300) 1.0 
=3.75 x 105 m 3 . 


Answer: The dimensions of the crater for a 300 
kt contact burst are 
R 
a 
' 
. 
. 


""'i07!= 1.2 WO·08 = 1.2 (300)0.08 = 1.89, 
a 
R. =:1.89 (Va) 1f3 =1.89 (3.75 X 105)1/3 
=137 meters, 
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and 
D 
. 


V. lt3 = 0.5 W-{).l2 = 0.5 (300)-{).l2 = 0.252, 
a 
D. = 0.252 (Va )1/3 = 0.252 (3.75 X 105 )1/3 


= 18.2 meters. 


Uncertainty. The dimensions of the apparent 
crater obtained for the above problem have a 
range of uncertainties that are dermed by the 
following: 
R 
1.1 WO.08 E:;; V. 1;3 E:;; 1.4 WO.08 


a 


125 EO; Ra E;;; 159 meters, 
and 


Da 
13 EO; V. 1/3 
EO; 25 meters. 
a 


The characteristics of the medium represent a 
major uncertainty in the crater data presented in 
Figure 2·743 through 2-75e. Details concerning 
material properties and geologic structure are 
not usually 
known about any particular site, 
but, in many ~ases small changes in these media 
characteristics can cause large changes in crater 
dimensions. The largest changes usually occur in 
the presence of an intersecting water table or a 
layered medium. Therefore; the range of uncer· 
tainty for a generic geology such as "wet soil" 
can be quite large, especiaUy in the prediction 
of the apparent crater depth. 


i 
Related 
Material: 
See 
paragraphs 
2-46 


through 2·50; see also paragraphs 1·14 and 1-33. 


"-' . 
. ) 
'- 


Problem 2·23. 
Crater Dimensions for a Shallow Burierl Burst 


and 


gi.;,. 


.';" 


~~. 


I 
I 


~. 
I 
I 
I 
' 


Figures 2-74a through 2"75e include curves 
that indicate the apparent ,crater volume. as a 
function of HOB for I kt explosion in various 
homogeneous media. 
Scaling. 
For 
Yields, other 
than 
I 
kt, 
the 
crater volume scales as follows:, 


~ 
-=W3a 
Val 


where Val is the apparent crater volume for a I 
kt explosion, and v;. is the corresponding vol- 
ume for a yield of Wkt. The height of burst 
.scales as 


Example: 
Gil'en. A hypothetical 30 kt burst in wet soft 
rock at a depth of burial of 3 meters. 
Find. The apparent crater radius and depth. 
SoluTion: The value of a is obtained from 
.Figure 2-76a, which must be entered with the 
scaled DOB. 
The scaled DOB is 


DOB ~ 
3 
- 0 965 
/k 1/3 
. Wl/3 - (35)i73-' 
m 
t 
. 


From Figure 2-76a, 


a = 0.299. 


The scaled HOB is 
HOB-3 
. 
HOB). =JiiCI ='(30)0.299 = - 1.09 mlktl1 .'" 


From Figure 2-75d" the apparent crater volume 
for a, I kt explosion, at HOB = -1.09 m/ktl1 in 
wet soft rock, is, 


VaJ=2.20x 
1Q4 m3/kt 3a 
, 


and 


v;. = Val (W)311 = 2.20x 104 (30)0.897 


= 4.65 x lOS m3. 


The' dirriensions for the 30 kt explosion are de- 
termined from the following expressions: 
. 
" 
R 
. 
-vm = 1.2 Wb. 


I 


'D 
.v.1iJ =05 W-C. 
• 


From, Figures 2-76b 'and 2-76c the values of the 
yield transition exponents, b and c, at a scaled 
DOB of I .09 m/ktl1 (HOB =-] :09 m/ktl1 ) are; 


b =0.023, 


c == 0.033. 


Answer: The dimensions of the crater for a 30 
kt explosion at a DOB of 3 meters are : 
R 
"'V'Th = 1.2 WO.023 = 1.2 (30)0 023 = 1.30, 
a 
Ra =1.30 (4.65 x 10S)I/3 '7,101 meters, 


and 


D 


Val,~3 = 0.5 W-O.033 = 0.5 (30tO.0 33 = 0.447, 


Da = 0.447(4.65 x '] OS )1/3 =' 34.6 meters. 


Uncertainty. 
The ranges of uncertainty for 
the radius and the'depth are 


] I WO.on lll;; JiL. EO; ] 4 WO.OB 
. 
Val/3 
" 


92 <;,R a <; 117 meters, 
'., 
D· 
. 


0.35 W..<l.033 <; y'1;3 <; 0.7 W- O.033, 
·a 


24 C;;D. <; 49 meters. 


Reklled Materi41: See paragraphs 246 through 
2-50; see alsoparBgraphs 1-14 and ]..33. 


-This could be written' DOB 1'"'!!2! .. 
3 
.. 1.09m/lul1 
. 
....,4 
(30)0.299 
' 
i.e., HOB"' -1.09 is eqlliYllall 10 DOB "'1.09. 
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\ , 


Problem 2-24. 
Crater Dimensions for a Low Yield Near-Surface Explosion 


and 


, . 


Figur~s 2-74a through 2,75e include curves 
that indicate the apparent crater volume asa 
function' of HOB for I· kt explosion in various 
homogeneous media. 


Sca1iJl8. 
For 
yields 
other 
than 
I 
kt, 
the 
crater volume scales as follows 


VI 
-=W3a 
Val 
. , 


where 
~l is the apparent crater volume for a 1 
kt explosion, and v.. is the corresponding vol- 
ume for a yield of k' kt. The height of burst. 
scales as· 


Example: 
Gil'en: A hypothetkal 3 kt burst over wet 
soil at a height of burst of I meter. 
Find: The apparent crater radius and depth. 
Solution: From paragraph 2-49, 


HOB 
= HOB = ..LQ... = 0.693 m/ktl/3. 
I 
I 
k'a 
(3)113 
, 


The apparent crater volume for a 3 kt burst will 
be determined by interpolating between the best 
estimate values from Figures 2-74b and 2-75b. 


~ (3,) = 0.6 Va (W <; 1) + 0.4 ~ (W > 1OJ. 


:-.;"... :-;0"';1'; :-74b, the apparent crater volume 
for a I. -kt explosion (for W <;; I), at a scaled 
HOB = 0.693 m/kt1/3, in wet soil is 


Val = 7.00 X 103 m3/kt3a , 
and 
VI (W":; 1),= Val (W)3a = 
7.00 x 103 (3)1.0 


. 
= 2.IOx 1()4 m3. 


FromF)gure 2-75b, the apparent crater volume 
for a 1 kt explosion (W > 10), at a scaled HOB' 
=0.693 m/k tl13• in wet soil is 
Val = I.S0x 103 m3/kt3a 


and 


V. (W> 10)=V.1 (W)3a= 1.80 x 103 (3)1.0 
= 5.40 x 103 m3. 
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Therefore, for a 3 kt burst 


V. =0.6(2:10x 1()4) +0.4 (5.40 x 103) 


= 1.48 x 1()4 m 3. 


Answer: The apparent crater dimensions for 
a 3 kt explosion with HOB of 1 meter are: 


Ra 
vm = 1.2 WO·08 = 1.2 (3)0.08 = 1.31, 
a 
" 
R a = 1.31 (Va )!/3 = 32.2 meters, 


, Da 
0.35 W-0.12 <; v. li3 <; 0.7 W-O,12 
• 
8 '" Da <;; 15 meters. 


If the 3 kt weapon in the preceding example 
were known to have a high ndiative output, 
similar to weapons of larger yield (10 kt), then 
the best estimate for the apparent crater volume 
would be found using only Figure 2-75b. There- 
fore, V. would be given by . 


Va = 5.40 X 103 m3, 


and 
'R. = 1.31 (V.)1/3 = 23.0 meters, 


D~ = 0.438 (V.)l/3 = 7.7 meters 
Related Materill1: See paragraphs 2-46 through 
2·50; see also paragraphs )-)4 and )-33. 


£, 
',,-,-' 


{ 


SUMMARY OF PROCEDURES FOR CALCULATING CRATERS FROM 
DEEP-BURIED BURSTS (DOBIW113 > 5 m/kt 113) 


... 
" J C"j 


"~ 


i 


f' 
II 
';. 
t" 


( 


Figures 
2-77 
through 
2-81,· together with 
appropriate scaling laws may be use.d to obtain 
crater volumes from deep-buried bursts in a vari- 
ety of situations. The following summary pro- 
vides a step-by-step procedure for obtaining such 
information. An example problem is provided 
immediately following the figures. 


Required (either given or estimated): 
Yield (W) in kt 
ActualDOB 


. Soil Type (one of the five generic types or a 
combination; (see Figures 2-82 and 2-86b). 


Yield-scaling parameters, a, for burst position: 
DOB/W lI3 > 5 m/kt 1/ 3 .... a = 1/3.4 


I. Compute DOB I 


(DOB for 1 kt) =Actu~/QDOB 


2. Determine 
~l (apparent-crater volume for 
1 kI),· using curve for given soil type from 
Figures 2-77 to 2-81. 


3. Compute 
Va 
(a"ctual· crater volume 
for 
yield other than 1 kt): 


Va (actual volume) = Val (W 3Q). 


Three curves are provided for each soil 
type: 
best 
estimate, 
upper 
bound, 
and 
lower bound. Using the foHowing guide- 
lines, select the curve most appropriate to 
the problem. 


a. Use best estimate v,J ;f no specific geo- 
logic data are given. If actual HE test are 
available for specific site, compare .them 
with the HE curves in Figure 
2-~2, and 
move toward the upper or lowe{bound 
NE curve to adjust for the.se data, de- 
pending on where the data are grouped. 


b. Use lower-bound Val' (if desirable) for 
target-oriented 
calculations 
(offensive- 
conservative) and to calculate Ra or Da . 
Consider Iower-bound 
Val for targeting 
when 
calculating 
crater-volume-related 
phenomena 
such 
as. ejecta, 
transient 


velocity~or displacement. ~ 


c. Use upper-bound 
Val (if desirable) for 
design-oriented 
calculations 
(defense- 
conservative) and to calculate R a or Da" 
.Consider 
upper-bound 
Val 
for design 
when calculating crater-related phenom- 
ena affected by volume. 


4. Calculate R a and Da using the following 
expressions: 


Ra =1.2 V.l/3 


Da =0.5 V.lf3 


The dimensions R a • Da• and V. can. be re- 
lated to other pertinent crater dimensions as 
follows: 


• 
The· radius to the crest of the apparent- 


c~ater lip is 


Ra~ =1.25 Ra· 


• 
The height of the apparent lip is 


0.20Da C; Ha~ <; 0.25 Da 
for deep-buried bursts. 


Variations in the thickness of the ejecta as a 
" function of range from the surface &founCi zero 
(SGZl are discussed in paragraphs 2-52 and 2-53. 


-In all instances it iJ recommended that IIPPer and lower bOllnd 
wall1es" be calCll1ated in order to lain an appreciation of the 
effecu of W1ccrtaiDtics on the particll1ar problem WIder con- 
lideration. 
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Figure 2-77. 
.trlter Volume IS • Function of Depth of Burill for I' Kiloton 
Explosion Buried in Dry Soil 
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Figure 2-78. 
Crlter Volume IS I Function of Depth of Bunl' for I , Kiloton 
Explosion Buried in Wet Soil 
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Figure 2-79. 
Crlter Volume IS I FunCtion of DePth of Bunll for I , Kiloton 
Explosion Buried in Dry Soft Rock 
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'Crlter Volume IS I Function of Depth of Burial for 11 Kiloton 
Explosion Buried in Wet Soft Rock 
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Figure 2-81. 
Crlter Volume IS I Function of Depth of Bunll for. 1 Kiloton 
Explosion Buried in H.rd Rock 
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Figure 2-82. 
Crlter Volume IS I Function of Depth of, Bunt ,for HE Explosions 
in Five Homogenous Geologies 
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Problem 2-25. 
Calcul,ation of Depth of Burial for Optimum Crater 


., 


Figures 2-77 through 2-81 show curves that· 
.indicate the apparent crater volume as a func- 
tion of DOB for a I kt explosion in various 
media which are indicated.in each figure. 


Scaling, For yields other than I kt the crater 
volume and depth of burial scale as 


DOB = WI/3.• 
DOB) 
• 
' 


.~ 
-= W3/3.4 
Val 


where 
DOB I and Val are the depth of burst and 
apparent crater volume for a I kl explosi<...I, and 
DOB and ~ are the, corresponding quantities for 
a yield of h' kt. 
Example. 


Gi~'en' A hypothetical 3 kt explosion in hard 
rock.. 
Find: The depth of burst at which the maxi- 
mum apparent crater will occur, andJhe dimen- 
sions of the' apparent crater at this depth of 
burst. 
Solution: From Figure 2-81, the maximum 
crater volume from a I kt explosion in hard 
rock occurs at a DOB I of 30.5 meters (note that 
the curve is flat in this -region, and a variation of 


~ few meters'in DOB· will have a small effect on' 
V). The apparent volume of the crater from such 
a burst is 


VI =7.20 X 104 m3, 


Answer: The corresponding depth of burst 


and crater volume for a 3 kt explosion are 


DOB.= (DOBl ) (W)1/3.4 = 30.5 (3)1/3 .• 


~ 42.1 meters, 


V. =(Val) (W)3/3.4 =7.20 x ] ()4 (3)3/3.4 
=1.90 x lOs m3 . 


.. The apparent Crater dimensions are defined by 
the expressions 
R. =1.2V.l/3 =1.2 (1.90 x 105 )1/3 
=69.0 meters, 


D. = 0.5 V.l/3 = 0.5 (1.90 x 105 )1/3 
=28.7 meters. 


Uncertainty,~The dimensions of the apparent 
crater have a range of uncertainty that is defined 
by the expressions 
1.1 val /3 <'R. <; 1.4 ~1/3, ' 


63 <;R a <; 80 meters; 


0.35 val /3 '" D. <; 0.70 V.l / 3, 


20'" D. C;;; 40 meters, 
In addition, the uncertainty in crater volume is 
given by 


V. (min) = Val (min) W3/3·4, 
V. (min) =2:80 x ]04 (3)3/3.4, 
= 7.38x 1()4 m3, 


V. (max) = Val (max) W3/3.4, 
V. (max) =1045 x 105 (3)3/34, 
=3.82 x lOS m 3. 


Related Material: See paragraph 2-5]; see also 
paragraph 1-14 and ]-33. 
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EJECTA 


. Crater !ejecta consists of soil or rock debris 
that is thrown beyond the boundaries of the. 
apparent crater by an explosion. Together with 
the fallback, which lies between the true and 
apparent crater boundaries, ejecta comprises all 
material that is completely dissociated from the 
parent medium by the explosion. The ejecta 
field is divided into two zones: 


• 
The crater lip, including the continuous 
ejecta surrounding the apparent crater. 


• 
The discontinuous ejecta, comprised of 


th~ 
discr.ete 
natural 
missiles 
that 
fall 


beyond the limit of continuous ejecta. 
The principal parameters used to describe the 
. ejecta are: 


• 
The average lip crest height. 


• 
Th'e radius to the crest of the crater lip 
from surface ground zero (SGZ). 


• 
The depth of deposition. 
, 


• 
The ejecta mass density. 


• 
The missile size distribution,which is a 
function of the azimuthal and radial dis- 
tance from SGZ, the radius to the outer 
boundary of continuous ejecta, and the 
maximum missile range. 


The principal variables that control the ejecta 
parameters are the yield and geometry, and the 
physical' characteristics of the earth medium. 
Figure 
2-83a 
shows 
the 
throwout {)f ejecta, 
from the SEDAN Event, a 100 kt explosion at a 
depth of 635 feet. 


2·52 
Ejecta Thickness 


The amount and linear extent of the continu- 
ously deposited ejecta in the crater lip are deter- 
mined primarily by the yield. The radial limit of 
continuous ejecta, which is the o~ter edge of the 
lip, will usually vary between 2 and 3 times the 
apparent crater radius. The lip crest height above 
the original surface is estimated to be one-fourth 
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of the apparent crater depth 
(Ha~ = 0.25 Da• 
paragraphs 2-50 and 2-51). The depth of ejecta 
decreases rapidly in an exponential fashion as 
the distance fromSGZ increases. About 80 to 
90 percent of the entire ejecta volume generally . 
is deposited within the area of continuous ejecta. 
Analysis of data for' craters formed by near· 
surface bursts over soil indicates that ejecta mass 
represents 
approximately 
55 
percent of the 
apparent crater mass (the remainder being found 
in fallback, compaction, and the dust cloud, 
which is blown away). When the same yields ex- 
plode at successively greater depths, ejecta mass 
increases significantly until the optimum DOB is 
reached. Ejecta thickness can be estimated for 


so~ as follo\"" (all dimensions are in feet): 


'e =O.9V. ~.6 R·3.86, for R>2.2V.1/3 


where ,e is th.e ejecta thick,,!ess, Va 
is the ap- 
parent crater volume, and Ris the distance from 
SGZ to the point of iitterest. This formula for 
ejecta thickness assumes that the ejecta mass 
density will be approximately equal to the ori- 
ginal in-situ density of the medium. The equa- 
tion may be considered valid for a soil medium' 
however, the bulking, which is inherent in a d~ 
turbed rock medium would result in greater 
ejecta thicknesses than predicted by the formula. 
Therefore, in a rock medium, the ejecli1 thick- 
ness should be mcreased by 30 percent to ac- 
count for the bulking, i.e., the formula for a 
rock medium is 


2-53 
Maximum Missile RInge 


Figure 
2-83b shows the maximum missile 
range as a function of depth of burst for a 1 kt 
explosion in two soil types. These data are based 
on empirical results from several nuclear ex- 
plosions. 
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Problem 2-26a. 
Calculation of Ejecta Thickness 


Paragraph 2-52 provides formulas for the 
prediction of ejecta thickness as a function of 
distance from surface ground zero and 'the 
apparen't crater volume (Va). Va may be ob· 
tained {rom Figures2.74a through 2·76c for' 
near surface bursts and 2·77 through 2·81 for 
deep-bu'ried bursts by the methods ,described 
in paragraphs 2-49 through 2·51. , 
' 


I 


F.xan:ple: 
Gh'en: A hypothetical 300 kt contact surface 
burst (center of gravity of the weapon 0.5 


, 
•• 
I 
meters above the surface) will produce a crater 
predicted to have an apparent radius of 137 
meters, an apparent depth of 18.2 meters and 
an apparent volume of 3.75 x lOS cubic meters 
in dry soil. (See problem 2-22). 
, 
' 
Find: The thickness of the ejecta at a distance 
of 300 meters from SGZ. 
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Solution: 


R =300 meters, 
'" 


2.2V 1/3 = 2.2(3.75x105)l/3 = 159 meters 
. 
. 
'R >2.2V.1/3 
Therefore the equation for soil given in para· 
graph 2-52 applies. 


Answer: The ejecta thickness is 
'e = 0.9V.1.6 R·3.86 , 
'e =0.9 (3.75xIOS)1.6 (300)-3.86 = 


0.20 meters 


Reliability: Based on empirical formulas de- 
rived from high explosive and nuclear burst. 
Specific reliability has not been estimated. 
Related Material: See paragraph 2-49 through 
2-52. See also Problem 2·22. 
. 
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Problem 2-26b. 
Calculation 'of Maximum Missile Range 
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, 
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Figure 2-83b shows the maximum range to 
which missiles might be expected from a I kt ex- 
plosion as a function of depth of burial in two 
soil types. 
Scaling For yields other than I kt, the maxi- 
mum 
missile range and depth of burst scale as 
fellows: 


where R mI is the maximum missile range for I 
kt, DOB, is the depth ofbuJ:St for I kt, and Rm 
and Dab are the corresponding range and burial 
depth for W kt. 
Example: 
Gil'en,' 
A hypothetical 3 kt explosion at a 
depth of 85 meters.in hard rock, 
Find' The maximum missile range for the 
explosion, 


Solution: The corresponding depth of burst 


for a I kt explosion is 
. 


DOBJ = DOB =~= 61 meters. 
WO.3 
(3)0.3 
From Figure 2-83b, the maximum missile range 
for a I kt explosion at a depth of 61 meters in 
,hard rock is 440 meters. 


Answer: The corresponding maximum missile 
'range fora 3 kt explosion at a depth of 85 
meters in hard rock is 


R m =R m1 WO.3 =(440) (3)°·3 =612 meters. 


Reliability: Based on empirical results from 
several nuclear explosions, but no specific relia- 
bility has been estimated. 


Related 
Material: 
See 
paragraphs 
246 


through 2-52. 
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Figure 2-83b. Maximum Miuile Range as a Function of Depth of Burial for a 
, kt Explosion 
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CHARGE STEMMING 


The 
tenn' charge 
stemming 
refers 
to 
the 
backfilling of material in the charge-emplacement 
hole. 
Ideally, 
charges should 
be 
completely 
stemmed and tamped to contain the explosive 
energy temporarily, which increaSes its coupling 
with the medium. Typical stemming materials 
are concrete, gravel, sand and water. The crater- 
ir.g curves shov.:n in preceding paragraphs pro- 
vide' predictions 
for 
fully 
stemmed 
charges. 
There are, however, operational considerations 
that may require reduced stemming or the capa- 
bility to emplace or remove the stemming ma- 
terial in a short period of time. Therefore, modi- 
fied 
stemming 
geometries, 
in~l~ding 
V'J;OUS 
degrees of stemming or no stemming at all, are 
importan t considerations in the evaluation of 
nuclear cratering phenomena. 
,The following changes in cratering phenom- 
ena generally may be expected when less than 
full stemming is used: 


• 
Air blast and the fraction of radioactive 
materials that vents will be increased. 


I 
• 
Energy loss out of the emplacement hole 
will reduce the coupling effectiveness of 
the explosion. This will result in crater 
dimensions that are smaller than those 
from fully stemmed explosions. 


• 
f::.nergy 
deposited 
in 
the 
emplacement 


hole 
will 
cause 
a modification to· the 
energy deposition pattern in the medium. 
The source will appear to be a distorted 
cylindrical source rather than' a concen· 
trated spherical source. Lip height, ejecta 
distribution, and maximum missile range 
will vary from those of a fully stemmed 
explosion. 


2·54 
Guidelines 
for 
Assessing 
the 
Effect 
of Stemming 


The 
crater 
dimensions 
from 
a 
partially 
stemmed or an unstemmed emplacement hole 


will be affected 
primarily by the amount of 
stemming, the 
diameter of the emplacemen t 
hole, and~the depth of burial. The following are 


some general conclusions concerning stemming 
that are based largely on HE experiments. 


• 
The crater radius, which is the dimension 
of greatest military importance, does not 
increase significantly with increased stem- 
ming. 


• 
Stemming about one-half of the emplace- 
ment hole (50% stemming) provides most 
of the crater depth that would be ex· 
pected from a fully stemmed charge. 


• 
Water 
appears 
to 
be a 'very 
efficient 
stemming material. 


Figures 2·84a, b, and c illustrate the HE experi- 
mental results concerning the effect of stem- 
ming on crater volume, diameter, and depth, 
respectively. 


.EFFECTS OF GEOLOGICAL FACTORS 


In most situations a weapon will not be burst 
in a homogeneous medium such as dry soil or wet 
soft rock. Instead, typical geologies may contain 
a water table at a shallow depth, a layering of 
one type of media over another (e.g., soil over 
rock), parallel planes of distant jointing in rock, 
or a steeply sloping ground surface. All of these 
factors can influence the formation of a crater 
and, 
in 
some 
cases, can change the size or 
characteristics of the crater significantly. The 
following paragraphs describe the general effects 
of geologic variations on cratering phenomena. 


2~55 
Sloping Topography 


Terrain slopes of about five degrees or more 
will affect the geometry of a crater fonned by 
either a surface or buried explosion. The ejecta 
distribution will also be affected. If the slope is 
gentle, the crater volume will be comparable 


2·193 


· -'., 
.:...:~- 
~'.::'.,'.~~..... :-.- 


3,.....-----.....--------,..------.....-----_ 


Z 
c SCALED 008, mlkt Y. 


w 
2 


~:lC 
...Jw 
O~ 
>~ 
a: w 
w~ 
~Ul 
l(z 
a:2 
UCw 
""::E 
O~ 
- 
Ow 
-~ 
, 


~Ul 
'- . 


a: 


100 
75 
50 


PERCENT STEMMING 


25 


ALL STEMMING MATERIALS 


OL... 
........ 
....L. 
,;,L,. 
......I 


o 


'k 


~"" 
I 


~ 
I: 
I 


Figure 2-84a. 
Increase in HE Crlter Volume IS I Function of Stemming for 
Vlrious DOBs 


2·194 
I 


3,..------"T"--------r-.....------.,------.., 


,- 
- 


Z • SCALED 008. mIlo Y. 


Z $2.5 
J 


~ 
! 
-Z < 3.5 


2.5 < Z $3.5 
1I, 


Cl: 
2!-'-------+---'-----+--"""":"-----i .----.--.. 
W 
~ 
W 0 
~ 
W 


C{ 
~ 
C ~ 
Cl: 
~ 
~~ 
C{ 
::J 
50 
W 


.... 
~ 
o 
~ 
OW 
~ to 
C{a: 
() 


. '1 


ALL STEMMING MATERIALS 


OL--------::~-----...I.----+-.&...o-----~ 
o 
25 
60 
75 
100 - 


PERCENT STEMMING 


Figure 2-84b. 
Increan in HE Crater Diameter IS a function of Stemmi"ll for 
Various DOBs 


2·195 


( 


Br----~I---~--~-----. 


Z • SCALED DOB. mllc1 y, 


2.5 )53.5 
, 


o..... 
A_LL_S-:T~EIM_M_~_G_M_A_T_E_R_IA_LS-:,:-I 
~i~{-_2_'6 
~ 


o 
25 
50 
75 
100 


PERCENT STEMMING 
'" 


z 
, 
0 
iiizc 
...... 
~~ 
-~ 
Q ... 


a:~ 
... 
1Il 
~Z 
«:l 


! £to 
Uw 
... 
~ 
O~ 
O~ 
~1Il 
, « 
a: 


," 


I 
Figure 2-84c. 
Increase in HE Crater Depth IS a FunC1ion of Stemming for 
Various DOBs 


2·196 


with that ofcraters on, level ground, but the re- 


, suiting crater will be asymmetrical; it will be 
"wider up-slopeand will have a larger lip down- 
slope, Figures 2-85a and b show, a photograph 
and a, sche.matic drawing, respe'ctively, of the 
crater from a nuclear event buried 30 meters 
below a 30 degree slope in tuff. The shape ,of the 
crater fonned by this explosion was influenced 
by the slope of the surface because almost' aU, 


, ' the debns fonned a rock slide originating at .the 


lower "Jge of the crater and terminating about 
240 meters down the slope. 


Limited 
small-scale' cratering' experiments 
have been conducted in moist, sandy soil and in 


de~ert alluvium on slopes ranging from 40 de- , 
grees to vertical wedges. For charges buried, on 
severe but nonvertical slopes, with DOB mea- 
sured normal to the sloping surface and with the 
verticai depth of overburden being greater thari, . 
.containment depth, crater dimensions decrease" 
with increasing slope, Optimum vertical DOB 
appears larger by about one-third than optimum 
DOB on level terraln, 


2-56 
Water Tables 


'A subsurface" ground water table 
i:h a soil 


medium will begin to influence the size and 
shape of the crater when its depth below the 
surface'is equal to or less ,than three-fourths the 


;.:~:::.:~cc apparent crater depth. Its effect is to 
flatten and widen the crater. Figure 2-86a illus- 
trates the effect of a shallow water table on the 
shape of the crater., fu the water 'table depth 
decreases, 
its 
effect 
becomes more evident. 
Based upon HE experiments, it can ,be con- 
,cluded that for a sunace detonation, the in- , 
fluence of the water table is not significant for 
scaled water table depths. greater, than about 
3.5 m/kt I/3 , and the apparent radius, volume 
and 
lip 
height increase exponentially as the 
water table approaches the ground surface. For 
shallow buried bursts (DOB/WI/) < 3 m/ktl/3 ) 


the fmal crater radius may be. as much as 50 per- 
cent greater than the predicted value and the 
depth as little as one~third ,that of the predicted 
value. , 


Figure 
2-86b shows a normalization tech- 
nique 'that has been developed for' plotting 
crater volume for surface and shallow 
buried~ 


explosions 
in a layered geology. 
A shallow 
.water table effectively causes a homogeneous 
medium to act as a two layered medium 


, d. == depth to water .table (thickness of upper 
layer) 


Va";,apparent crater volume 


~ =apparent crater volume predicted in the 


upper layer alone (usually dry soil) 


VL =apparent. crater volume predicted in the 


lower layer'alone (usually wet soii), 


The curve fitted t~ the water table crater volume 
data in,Figure 2-86b is givenby the expression 


v.-v.'" 
' 
V. l =l-exp(-S.4dlV.1(3). 
u - L, 
' 


- , 
' . 
Iteration is required to solve this expression for 
v.; 
however, 
the 
technique 
converges 
very 
rapidly: It is recommended ,that' the curve in 
Figure 2-86b be used to predict crater volume 
for sites having a shallow water table. 


2-57, Bedrock 


The influence of a bedrock layer below a soil 
,medium is similar to that of a, water table, 
though somewhat less pronounced. For surface 
explosions it appears that a hard layer may limit 


, the crater depth wHile causing the crater shape 


to be more nearly trapezoidal than parabolic or 
conical (see Figure 2-86c). The bedrock layer 
may increase the crater radius slightly (5 to 10 
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Figure 2-851. 
NEPTUNE CraUr and its PClIition on Slope of MISa 
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Figure 2-86b. 
Crllering0.,. for Dry Soil over WeI Soil 
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percent).:The layer may decrease the fmal crater 
depth bYl as much as one-third when the over- 
burden layer is as shallow as one-fourth the pre- 
dicted apparent crater depth. 


i 


Although the effect of a water table and the 
effect of a hard layer are different, the water 
table curve is a good approximation for predict- 
ing the crater 
voh.~me in a layered geology in- 
volving a gradual transition ·to rock. The water 
table curve should not be used for geologies 
involving a sharp transItion to rock. 


, 
The two-layer problem can be generalized to 
three layers. Let VI' V2, and V3 be the crater vol· 
ume associated with.the top, ITiiddle, arid bottom 
materials considered separately. Then, V23 ii the 
crater volum~ associated with tne middle layer 
and bottom, layer, and V is the crater volume in 
the total layered medium. 
V23 is calculated by neglecting the top layer 
and then solving a two-layer problem with V2 as 
the upper layer and V3 
a~ the lower layer. The 


thickness of the middle. layer d 2 is used as the 
depth t~ the lower layer. This ~ves 
. 
, 
V23 - V3 
. 
= l-exP(-5.4d2/~31/3) 
.~ - V3 


" 
;,. 
t ....... c-aJculated by solving another two- 
layer problem with VI as tlle upper layer and V23 
as the lower layer. Tne thickness of the upper 
layer d I is used as the depth of the lower layer. 
This gives 


. V- V23 _. 
3 


V, 
V 
- l-exp(-5.4d l /VU ). 
1- 
23 


2-58 
Rock Beddina!Jointing 


If a low-yield explosion occurs at or very near 
the surface, the bedding or joh,ting planes in 


Tt)C~ ra.'! influence the shape of the crater pro:- 
duced. The direction of the ejection process 
I 
2-202 


will also be affected. The formation of the crater 
will tend to follow the direction of the pre- 
dominant joints, and the crater radius will be 
increased by as much IS one-third in the direc· 
tion parallel to the joints and decreased by as 
much as one-third in the direction normal to the 
joints. The magnitude of the crater depth is 
usually not affected significantly, but the deep- 
est point may be shifted to one side of the 
crater. As the yield or the DOB is increased, the 
influence of rock jointing is reduced. 
The dip of bedding planes will influence 
energy propagation, and this will cause the maxi- 
mum crater depth to be offset in the down-dip 
direction. Little overall effect is noted on the 
size of the crater radius, but differences in ejec- 
tion angles cause the maximum lip height and 
ejecta radius to occur down-dip. 
. 


2-59 
Snow and.lce 


Measured craters in snow or ice are a rarity; 
however, data for a few craters have been re- 
corded for surface HE explosions. In general, 
these 
craters are 
larger than would be pre- 
dicted in soil and they are characteristically wide 
and flat. The trends in crater size and shape for 
a surface explosion in a snow/ice medium are 
given by the expressions 
R. =38.8 W0.26 meters 


D. =5.75 WO.IS meters 


where W is the yield in kilotons: 


MULTIPLE BURST GEOMETRIES 


Nuclear weapons may be detonated in close 
proximity to each other to create a linear crater 
or a series of interconnected craters. When 
simultaneous explosions are sufficiently close to 
cause interactions, the shot geometry is termed 
umultiple burst." A linear array of this nature is 
referred to as a 
lOW crater, as illustrated in 
Figure 2-87. 
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Figure 2-8&. 
EHecti of a Hard Near·Surf~ LIVer on em-ring 
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Figure 2·87. 
A Row Crater Produced by Simultaneous 
Exr'osions of F:ve 1.1 kt Nuclear Devices 
at Optimum DaB (BUGGY Event) in 
Basalt; DOB= 41 m (135 ft.); Spacing = 
45.7 m (150 ft.); Width =76.2 m 
(250 ft.); Depth =19.8 m (65 ft.); 
Length = 262 m (860 ft.) 


In 
a~dition to the parameterS normally asso- 
ciated with an underground explosion, the re- 
sults of a row charge also depend upon the 


spacing between charges (s) and the degree of 
simultaneity of the explosions. Differences of 
only a few milliseconds in the times of explosion 
of adjacent charges can result in significantly de- 


lOl'CUCU ruw charge crater dimensions. Close spac- 
ing of weapons (less than 1.4 times a single 
weapon crater radius) increases both the crater 
radius and depth compared to single weapon, cra- 
ters at the same scaled DOB. Experience has 
shown ptat a spacing of 1.0 R. (for explosions 
atoptimurn DOB) results in a smooth channel 
with 
cr~ter depth and radius approximately 20 
percent greater than those of single explosions. 
A spac~g of 1.25 R. will still form a smooth 
channel, but with less enhancement of crater di- 
mensions, while spacing of 1.5 R. will provide 
an adequate linear obstacle. The length of a row 
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crater L can be found in the equation; 


L=s(n-I)+2R., 


where n is the number of charges in the row, and 
s is the spacing between charges. 


SECTION III 
II 
GROUND SHOCK PHENOMENA. 


2·60 
General Nature of Free-Field 
Effects 
_ 


• 
The position of the weapon' at the time 
of 
detona~i')n has ;, marked influence on the 
amount and nature of the energy transfer into 
the ground. The effects at depths below the sur- 
face from a high air burst will be primarily of 
the air-induced type. For surface or near surface 
bursts, both air-induced and direct-transmitted 
ground shoc:k effects are produced, and both 
may be of importance. Deep underground bursts 
cause dire c t-u an smitt e d effects almost ex- 
clusively. 


• 
The nuclear field tests of the past several 
decades involving air, surface, and fUlly contain- 
ed bursts have provided the data that permit 
development of the expressions for predicting 
the 
peak 
values 
of 
air-induced 
and 
direct 
-transmitted acceleration, velocity, and displace- 
ment. Field observations for regions where the 
air blast overpressures are less than 1,000 psi 
generally indicate that, within about 50 to 100 
feet of the surface, ground motions are predomi- 
nantly air-induced. The relatively few measure- 
ments of acceleration and displacement that 
have been made at significant distances below 
the ground surface show the effects of attenua- 
tion 
and 
dispersion 
as the air-induced pulse 
travels through the ground. 


_ 
A typical acceleration-time record from 
a gage mounted near the surface of the ground 


_. 
. 
\ 
~I 


_ 


• 


Reproduced from 
best available COpy. 


'---~---- 


• 
exhihits a systematic pulse shape corresponding 
to the passage of the shock wave, with a random 
type 
disturbance often superimposed 
on 
this 
systematic pulse. In the superseismic case, when 
the air blast propagation velocity is greater than 
the ground seismic velocity , the downward ac- 
celeration is large compared with the following 
upward acceleration. The velocity-time record is 
similar in shape to the overpressure-time pulse, 
at least in the early stages. In cases in which the 
ground motion outruns the air blast, there is a 
slower rise in the acceleration; it may even be 
reversed In direction inItially, i.e., upward rather 
than downward. In such cases, the signals may 
last for a longer time than the positive phase 
overpressure 
duration. The ussociated velocity 
wave form usually exhibits a velocity jump as 
the air-shock wave passes over the position. but 
the overall record is characterized by a consider- 
able degree of oscillation. When precursors are 
present,the acceleration and velocity records 
may exhibit higher frequency components and 
more random type of oscillation. Wave front dia- 


Compression 


WOVI 


(A) Superseismic 
Blast 
Wave 


grams for the superseismic and outrunning cases 
are shown in Figure 2-88. 
. 


• 
Most values of velocity. have been ob- 


ta= by 'integrating acceleration records. In 
general, velocities appear to be more predictable 
quantities than accelerations or displacements. 
Displacement data are even scarcer than velocity 
data, and only a few direct measurements have 
been made. 
For the most part, displacement 
data consist of integrated acceleration records or 
extrapolation of displacement spectra to zero 
freJenCy. 
' 
. 
Knowledge of direct-transmitted ground 
shoc 
effects is limited. Those data that exist 
were obtained from buried high explosive shots,. 
from 
several surface and near surface nuclear' 
detonations, and from a number of fully con- 
tained 
bursts. For .a completely buried shot, 
measurements indicate that at any particular 
range R. the strain or velocity. pulse rises in a 
parabolic manner to a peak value in a time 
roughly corresponding to Tr =R/6cp to R/12cp , 
where R is the range from the explosion, and cp 


(B) 
OutruMing 
Ground 
Wove 


Figure 2-88.. 11 Wave Front Diagram for SuiiismiC Air 
Blast and Outrunning Ground Wave 
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II 
is the effective seismic velocity of the medium. 
The signal then drops slowly with a total posi- 
tive phase duration corresponding approximate- 
ly to To =R/2cp to Rlcp ' The effects of layer- 
ing, reflection', and refraction, tend to introduce 
high frequency ,and random components of mo- 
tion. 


2-61 
Types of Effects 
_ 
Ira it is convenient to' consider the earth 
shock resulting from a nuclear explosion as pro- 
ducing both. systematic and 
random effects. 
Systematic effects can be. divided further into· 
two major ty~es: 


• 
Air-induced 
shock 
associated 
with 
the 
passage of an air shock wave -over the sur- 
face of the ground, and the overpressure at 
the surface above the structure transmitted 
downward with such attenuation and dis- 
persion 'as may be consistent with the 
physical' conditions at the site. 
• 
Direct-transmitted ground shock that arises 
from duect energy transfer' from surface, 
near surface, or underground bursts. 


_ 
Random effects include high frequency 
ground-transmitted shock, surface-wave effects, 
reflections, refractions, etc. The dominant effect 
depends on such facto~ as weapon yield, point 
of detonation with respect' to the ground sur- . 
lace, range from ground zero, depth of the mea- 
surement, 
and, 
in 
particular, 
the 
geologic 
conditions. 
i 
tit Reasonable estimate~ can be made of the 
maXlmum values of displacement, velocity, and 
acceleration 
associated 
with 
the 
air-induced 
shock nnd in more restricted cases for the direct- 
transmitted ground shock under more or less 
uniform geologic conditions. 


2·62 
Air·lncsuc:ed Effects II 


• 
At the surface in a homogeneous medi- 


um~e maximum transient elastic vertical dis- 
placement may be expressed as 


Hpl/) 
drc = 2E""' 


where H is the depth to which the air-induced 
shock wave extends during the effective dura- 
tion of the shock, t.. P 
is the side-on overpres- 
I 
10 
sure at the surface. and E is the restrained mod- 
ulus. The dUration of the shock may be approxi- 
mated by 


for 


2 <; P10 .a;;;) 0',000 psi. 


Thus, 


The restrained modulus is related to the effec- 
tive 
seismic velocity, 
the 
mass density, and 
Young's modulus of elasticity of the medium as 
follows: 


E = c2 =(E) l 
I - 
11 
) 
P 
p 
p \< I + 11)(1 - 
211) , 


where E is the Young's modulus of elasticity 
(psi), p is· the mass density of the medium (P = 
..,Ig Ib sec2Jft4 , where "I is the unit weight of the 
medium in pounds per cubic foot and g is the 
gravitational acceleration 
constant), and 
1/ is 
Poisson's ratio. For values of 1/ of 0.25 or less, 
the relation may be approximated by: 


E = PC~ - 
E. 


If the unit weight of the soil is about 115 Ib/ft] • 


~ 
.- 
an approximate value of E is 


E = 25,000 C,~~oy, psi.· 


Substituting the values for H andE gives 


II The permanent vertical, displacement i.s 
not so easily estimated; at present the best estI- 
mates can 
be made by considering the static 
stress-strain properties of the soil. In the absence 
of better information. it is suggested that the 
permanent displacement at the surface be taken 
to be· the 
following 
for overpressures greater 
than 40 psi: 


_ Pso 
- 
40 ( I ",000 )2 . 


dsp 
- 
.30 
-c- 
mches,. 
p 


where c 
is the seismic velocity of the material 
p 
. 
near the surface. 


• 
From wave propagation theory it may 
be ~onstrated that the following relationship 
exists for particle velocities in terms of. the di- 
. mensional units defined above. 


where 
£ 
is the strain. Upon substitution this 
gives the maXlmum v.elociiy 


( Pso) (I ,000) .. / 
vs = 50 
100 ---c;- 
m. sec. 


_ 
For a rise 
time of velocity of about 
0.002 seconds, and considering a parabolic rise 


with a 20 perceht,increase for nonlinearity, the 
maximum acceleration at the surface. computed 
as a rate of change of velocity, may be expressed 
by 


( Pro)(I ,000) 
,as 
:::= ISO 
100 
~,g. 


_ 
Since the surface acceleration is not re- 
lated. solely to the maximum velocity, but is par- 
tially dependent. on the weapon yield and other 
factors,. it is recommended that a value of cp no 
greater than 2,000 fps be used in this expression. 
re=ess of the actual surface seismic velocity. 
_ 
Horizontal effects data are lacking for 
the surface. At pnisen~, it is recommended that 
the maximum horizontal deflection be taken as 
one-third the vertical, the maximum horizontal 
velocity as two-thirds the vertical,and the maxi- 
mum horizontal acceleration as equa:I to the ver- 
tical." 
. 
. 
• 
The same type of reasoning can be fol- 


Io: in arriving at maximum values of accelera- . 
tion, velocity, and displacement at depths below 
the surface. The difference in deflection be- 
tween the surface and some point at a depth y 
beneath the surface cannot exceed the surface 
.stress divided by the lower restrained modulus 
of deformation, E, in the interval y, and multi- 
plied by the depth y. Thus between the surface 
and a depth Y. not greater than 100 ft, assuming 
no attenuation of pressure with depth, an upper 
limit of the elastic component of differential dis- 
placement is given by 


_ 
The actual difference in denection· i 
more likely to be one-half this value, and may b(; 


If the unit weight of the medium differs significantly from 
1 
Ib/ft3• a reprelentative value of'l should be used rather than 
the approximation liven here. 
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I 
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II 
comidered to vary linearly to a depth of 100 ft, 
which results in the following expression 


d. - dye = 2.4(;06)~1 ,~or C~·rt)inCheS 


I, 
P, 


I 
,where y. = y for depths of 0 to 100ft, and y. = ' 
100 ft for depths greater than 100ft. 
II It is r:t0t generally considered desirable 
to assume any attenuation of displacement, ex- 
cept for very soft soils. Sirice total impulse must 
be preserved, the wavelength ,of the strain pulse 
probably. incr~ases with depth, and, therefore, it 
does not seem ,reasonable that there shOUld be a 
reduction in elastic component of displacement 
consistent with the reduction in peak stress, as 


--given by the 
~ttenuation relationships for stress 


orldty with depth. 
' 
, 
It is recommended that the permanent 


disp acement of soils at depths below the surface 
be 
dete~ed 
by multiplying the permanent 
displacement at the surface by (lOa - y·)/IOO; 
i.e., 


, 
( 
, 
d 
= 
100 - y. 
YP 
d~ 
!IOO) 


_ PIC) - 40 II,OOOY(lOO-y.);, 
' 
--"3"o-\-C;-J 
100 
mches. 


Thus, permanent displacements belowJOO feet 
would be taken to be zero. At depths less than 
100 feet, the total displacement should be taken 
as the sum of the elastic displacement at the 
surface, die' shown previously, and the perma- 
nent displacement, d p' given above. 
• 
The value ormaximum stress below the 


~, is given approxirriately, for any yield or 
overpressure; by the following expression 


Py = G'p1O 


, 2-'15 
Z#. 


where 


, 
1. 
CIt =---. 


, 
I +-f 


and. for depths that are not near the surface 


"H = 230\ (100\0.5 (W(Mt)) 1/3. 


PIO J ' 


_ 
Note: The attenuation factor, a, that is 
used in this section is different than the corre- 
sponding attenuation factor 
l/~ that is used in 
Section II. "Damage 
to 
Belowground Struc- 


tures," of Chapter 11, "Damage to Structures." 
The. attenuation factor given in Figure 11-29, 
Sechon II of Chapter 11 is the result of more 


recent studies. It reflects the influence of soil 
type (as represented by its seismic velocity) on 
peak pressure or soil particle velocity attenua- 
tion. The disparity between the two factors in- 
creases as the peak ove,rpressure and/or the soil 
seismic velocity increase, but the disparity is rel- 
atively insensitive to variations in weapon yield. 
As a result of the influence,of seismic velocity, 
the differenCe between the two factors also in- 
creases as the depth below ground surface in- 
creases. The factor given in Figure 11-29 gives 
substantially 
smaller attenuation rates under 


high overpressures and high seismic velocities 
than does the factor used in this section. Conse- 
quently, from the "vulnerability analysis" point 
of view, the factor of Figure 11-29 is more con· 
servative, and its use in shock analysis studies is 
generally recommended. The attenuation factor 
obtained from Figure 11-29 may be introduced 
into the equations of this section by simply re- 
placing "a" of this section with "lIP" where 
"P" is obtained from Figure 11-29. As ~ention­ 
ed previously, it is recommended that this re- 
placement be made, in aeneral; however, there is 


where t r is the effective velocity pulse rise time. ' 
This expression gives values which appear to 
agree well with test data for depths greater than 
lot; 
, 
- 
In the absence of better information at 
dep 
s below the surface, it is recommended 
that the ratios of horizontal to vertical effects be 
taken equal to those previously specified for the 
sUrface, Le., horizontal deflection is one-third 


v y = Ct:Vs' 


_ 
The peak downward acceleration attenu~ 


ates sharply with depth, especially near the sur- 
face. The sharp attenuation is caused primarily 
by the increase of rise time in the stress or veloc- 
ity with depth. 
, 
e 
The time of rise of the maximum veloc~ . 


ity from an initial zero value can be taken as 
, one-half the transit time of the shock wave from 
_the surface to depth considered. For a parabolic 
rise in velocity the acceleration can be consider- 
ed to be twice the value for a linear rise, which 
"leads to the expression - 


-' 
one important exce,ption. The factor 1/{3 is not 
compatible 
with 
the 
shortcut 
vulnerability 
curves shown ,in Figures 11-48 thrcugh 11-50, 
and calculations of similar curves that are com- 
patible with the 1//3 attenuation factor have not 
been performed at this time. It is for this reason 
that the ex attenuation factor is introduced here, 
and its use will be illustrated ,in Problem 2-27. 
The use of the I f(3 attenuation factor is illus- 
trated in Problems 11-3 and 11-4. 
I 
At depths below the surface, it is,recom- 
men ed that velocities be attenuated in about 
the saine manner as the maximum stress. Thus, 
the vertical velocity at any depth y is given by 


P - 
2 
p 
t 
- 
1,+ 
Vi 
i 


the vertical, the maximum horizontal velocity is 
two-thirds the vertical, and the maximum hori- 
zontal acceleration is equal to the vertical. 
_ 
Layered media pose a complicated situa- 


tio~ut with care and judgment, reasonable 


.' estimates ofaceeleration, velocity', and displace- 
ment can be made. A convenient method is to 
use the stepwise passage of a stress wave down- 
ward through the medium. In general, the basic 
concepts governing the computation of displace- 
ment, velocity. and acceleration are the same as 
those described for a uniform medium. The dis- 
placements at any particular time may be com- 
puted by dividing the average preSsure in an in- 
terval by the modulus of elasticity to obtain the 
strain and then multiplying by the length of the 
interval to obtain the displacement; the total dis- 
placement occurring over the length of the pulse 
is ~um of the incremental displacements. 
' 
• 
Complications arise at the interface of 
two media because of stress transmission and 
reflection. For soil and rock, the interface may 
not be sharply defmed and the reflected and 
transmitted stresses probably do not follow the 
'laws governing purely elastic media. 
If it is 
known that the interface is fairly sharp, an esti- 
mate of the reflected and transmitted stresses 
can be made from the following relationships: 


where tV is 'the ratio of the impedances of. the 
two media. The stresses at the interface must be 
taken as equal, and, from considerations of con- 
tinuity, the displacements also are equal. 


,2-63 
Outrunning Ground Motion' II ' 


_ 
Outrunning ground motion occurs when 
the air-shock velocity U decays ,below the sei$. 


= 5 (t~~) ~ ex, g, 
Y 
'. 


v 
= 2..!. 
0y 
't 


r 


_ 
Without 
seri~us error, th~ relation. be- 


twPoP" oealntrain e and peak partIcle veloCIty v 
is give~ by the equation previously shown in 


II 


mic wave velocity c of the medium.In the most 
general 
s~nse, outrunning ground motions and 
direct-transmitted 
ground shock are different 
phenomena. ,Direct-transmitted ground shock is' 
ground motion propagated through the ground 
media from the region of the crater (this is often 
denoted outrunning when it clearly arrives, as 
reflected or' refracted waves from deep-seated 


'~Vf"r<;, 
~he::ld of the blast front). Outrunning 
gr'ound' motion 
occurs when 
the air-induced 


, ground motion begins to propagate more rapidly 
than the air. blast shock. front, as illustrated in 
Figure 2-88, In actual fact, the motion-time his- 
tory at a pointin the medium 'can be quite com- 
plicate:', 
even 
random in 
nature, when 
air- 
induced, outrunning, and direct-transmitted ef- 
feil.i.llI arrive at about the same time,' 
,- 
_..' • 
Some infonnation on outrunning' m~ 


tions observed in field tests can be found m 
"Nuclear Geoplosics" (see bibliography). Com- 
,puter codes that are presently being developed 
should provide guides to even better estimates or 
outrunning motions. 


2-64 
Direct-Transmitted Ground Shock tI 


_ 
The .energy transmitted directly to the 
earth from a surface or near-surface burst can be 
propagated effectively through competent ma- 
terial for long distances. Experimental data that 


~......umu .. i.~ .his are available from buried high 
explosive shots and from surface and contained 
nuclear weapon tests. 
_ 
Fora completely buried shot, the fITSt, 
po=' of the strain o'r velocity record at a dis- 
tance R from the point of burst has the form 
described in paragraph 2-60, i.e., 


o 


2vr 
12c~ 
24c~' 


Dr =-- ==: -- f to"-R" f. 
Tr 
R 


paragraph 2-62, i.e., 


. 2 
I' 
2 
dr =~vrTo ==:T Re tOTRe. 


_ 
If data for strain,acceleration, or dis- 
placement are available, approximate relations 
for, the other quantities may be obtained. In gen- 
eral, the seismic velocity enters into the relation-" 
shil' 
s as indicated in equations shown above. 
Most of the available test data from 
whic 
direct-transmitted shock effects may be 
estimated were obtained from buried nuclear 
and high explosive detonations. To extrapolate 
from the test data toestirnate shock effects pro- 
duced by surface nuclear detonations requires 
establishment of the equivalence factor relating 
buried HE to buried nuclear yield, and also 
establishment of the equivalence of fully buried 
detonations to surface detonations of the same 


ty~fexplosive. 
' 
A yield effectiveness factor of 0.2 to O.S 
has 
n used to' relate fully contained nuclear 


, _ 
An estimate of the, maximum displ.ace- 
ment d,.jn the radial direction may 
b~ .obtamed 
by integrating the ·area under the .posltlve phase 
of the velocity curve. If the velocity waveform is 
parabolic, then 


"Assuming that the steepest part of the velocity- 
time curve has a slope that is twice the average 
slope during the rise phase leads to the following 
'relationship between the range of peak accelera- 
tion, peak velocity, and strain, all in the radial 
direction. 


R 
'R 
T· =--to-. 
o 
2cp 
cp 
and 
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. explosions to fully contained HE explosions. To 
relate a surface nuclear burst to a fully contain- 
ed nuclear burst, factors of 0.02 for a burst 
slightly above the surface to 0.05 for a burst 
slightly below the surface are recommended. Sig- 
nifican t 
uncertainties 
exist 
regarding these 
equivalence factors. 


_ 
Expressions 
for 
acceleration, velocity, 
and displacement have been derived on the basis 
of test data available for accelerations from shot 
RAINIER in Operation HARDTACK in volcanic 
tuff with a seismic velocity of 6,000 fps. For 
scaled ranges of 2.500 ft/Mt l /3 and closer, and 
an equivalence factor of 0.05 as just discussed, 
the acceleration may be' expressed as follows: 


Gr 
= 0.4 
(~~/(Mt))5/6 
(l.O~OY.5 C,~~oy, g. 


_ 
Using the relations between acceleration, 
velocity, and displacement'. and choosing the 
equations with cgefficients of 12 and I(3, re- 
spectively, this equation for acceleration may be 
used 
to arrive at the following equations for 
radial velocity and displac;ement for materials 
with seismic velocities less than 10,000 feet per 
second: 


_ I" IIW(M )\s/6 (1 ,000)2.5 ( 
cp 
). 
I 
vr 
- 
~ '\ 
t 'J 
~ .] ,000 
m. sec 


d~ =4 (W(Mt~S/6 (1,~0).5 inches. 
tit.There is Little or no information on tan- 
gential motions, and until additional informa- 
tion becomes available, it is recommended that 
maximum values of tangential acceleration, ve- 
locity, and displacement be taken as 1,2(3, and 
113 times the corresponding radial values. 


.lilt is desirable' to consider the direct- 
transmitted effects as being applicable only at 
some distance below grol.!nd surface, except pos- 
sibly at close in ranges. At scaled distances great- 
er than 
2,500 ft(Mt l / J , the limited available 
data suggest .a decrease in the rate of decay of 
. motion with distance. Approximate expressions 
for radial acceleration, velocity. and displace- 
ment for materials with an .average. seismic veloc- 
ity less than 10,000 fps are as follows: 


~ 
~ 
(,. 
) 75 ( 
C 
)2 
. 
_ 
0.25 
1,000 
. 
p 
. G r 
- 0.081. W(Mt) 
-R- 
1,000' g, 


II 
~ 
(1 000)1.5 ( 
C 
\ 
. 
vr = 4.8,W(Mt'l.0.s ~ 
1,0~0/n./sec 


( 
\ 
( 1000)1.25 
d r =3.19 
W(MOt75 R 
inches. 


These expressions give, only very approximate 
values. Geologic conditions can produce large 
and random variations in the motions at large 
distances. For materials with seismic velocities 
greater than 10,000 fps, e.g., granite rock rna· 
terials, the field test data suggest the following 
relationships as being applicable at scaled ranges 
of about 2,700 ft/Mt l / J and closer. 
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eLf) 


dr = 12(W(M't))S/6 
(1,~0t 
2 


~ 


165)5/6 (18,000)5/3. 
' 
-- 
. 
, mches. 


,oy 
cp 


II As would be expeCted from the different 


scaling relationships,. ~he expression for materials 
.,,'ith 
~~i~mic. vll'.locities less than 10,000 feet per 
second and those for materials with seismic ve- 
locities greater than 10,000 feet per second 
show lack of agreement at aseismic velocity of 
10,000 feet' per second, which divides the areas 
of applicability. This matter is currently receiv- 
i;;g further attention, but in the inter:J., it is 
believed that the expressions provide the best 
es.'tes or:such motions. 


-'. 
No expressions exist for motions at large 
distances in hard materials that arec9mparable 
to the expressions presented above for soft ma- 
terials, It was necessary to make approximate 
estimates to develop the relationships presented 
in Section 111, Chapter 11. 


2-65 
Spectrum Concepts - 
Simple Systems 
_ 


• 
Whl;n structural systems or equipment' 
are subject to a base disturbance, such as that 
arising from the ground motion associated with 
a nuclear explosion, the response of the system 
is governed by the distribution and magnitudes 
of the masses and resistance elements. A knowl- 
edge of the response of systems subjected to 
such loadings is extremely important from the 
standpoint :of design in order to proportion the 
.structure ~ that it will not undergo complete 
collapse, and to protect the structure, equip- 
ment, and personnel from shock damage. One of 
the simplest interpretations of ground motion 
data involves the concept of the response spec- 
trum, which is a plot of the maximum response 
of a simple linear oscillator subjected to a given 
input motion as a'function of frequency. 


2-"t84 
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_ 
The time history of the actual motions 
of the earth caused by the passage of a shock 
wave over the surface is complex and is subj~ct 
to considerable uncertainty. However, the prin- 
cipal effects on equipment and structural com- 
ponents can be described by use of the concept 
ofashock-response spectrum. 
. 
• 
Consider a piece of equipment or an in- 
ternal element of a structure supported at a 
point on an underground structure that is sub- 
jected to motion from blast. The equipment or 
element can be represented as a simple oscil- 
lator, consisting of a mass m attached by a linear 
spring of spring 'constant k to the base as shown ' 
in Figure 2-89(a). The natural frequency f of the 
oscillator is 


f =, w/2fr = _l_J k .- 


2fr 
m 


,_ For a given transient ground motion x, 
the mass m will be set into motion. If the mo- 
'tion of the mass m is designated by the coordi- 
nate y and the motion of the base by the coordi- 
nate x, the deformation in the ,spring u is (y-x) 
and is positive when y exceeds x. For such a 
system the most useful spectra are found to be 
the following: the maximum relative displace- 
ment u of the system; the maximum relative ve- 
locity u; the maximum pseudo relative velocity 
or the quantity wu. the circular frequency w 
-times the displacement u. which is not quite the 
same as uand differs considerably from uat low 
frequencies; and the absolute acceleration of the 
mass, w2U, which is nearly the same as y. A plot 
of the maximum values (denoted D. V. V. and Ag, 
respectively) as a function of the frequency of 
the system fot a given type of input constitutes 
a response spectrum for the particular function. 


, These quantities are dermed as follows: 


Ag: =Iw2DI = (2ft/)2D 


V =I c.JD1 
I: 2r[D, 
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Figure 2-89. 
_ 
Single-Oegree-of-Freedom System 
and T.ypical Response Spectra _ 


2~' 
VJ 


The maximum energy absorbed in the spring is 
associated with V. The quantity j7 is not particu- 


lar~seful. 


• 
For inputs of the' type associated with 
air-mduced or direct-transmitted ground shock, 
the spectra for a simple system will have the 
characteristic 
appearances 
shown 
in 
Figure 
2-89(b). The. limiting conditions that apply to 
the spectra illustrated in Figure· 2-89(b) are of 
particular interest. 
For example, as 
th~ fre- 
quency of the system approaches zero, the mass 
m does not move when the base of the system 
has a given 'motion applied. Then the relative 
displacement u is equal to the negative value of 
x. and the relative velocity u is equal to minus 
the maximum base velocity x. The dotted line irl 
Figure 2·89(b) indicates the difference between 
urn ax and V that may be observed in the region 
,'of low frequencies; likewise, there is some differ- 
ence between urn ax and Vat the high frequency 
end, but these differences are of a smaller order 
of magnitude. Typical values of V that might be 
expected for blast conditions are shown in Fig- 
ure 2·89(b). An additional control is that as the 
frequency of the system becomes large, the dis- 
placement y approaches the displacement x, i.e., 
u approaches zero, and the 
ac~leration j.imax 


approaches xrn ax' These controls are of special 
i~ance in arriving at design shock spectra. 
_ 
Ordinarily, t~e input for ground motion 
consIsts 
01 two parts, a systematic portion on 
which is superimposed a series of random oscilla- 
tions. The magnitude of the peaks of the ran· 
dom components may be either small or large 
compared to the systematic portion. The ran- 
dom part may exist over the entire range of the 
systematic portion, over only part of the range, 
or even may be prior to the systematic portion. 


_ 
For a random series of pulses, the rela- 
tive velocity peak of the spectrum compared 
with the maximum input velocity can be high, 
but it is not likely to be much higher than about 


~, '..nlt:S3 an almost resonant condition is obtain- 


ed with several pulses of alternate positive and 
negative signs of exactly the same shape and du- 
ration. Stich a resonant condition for velocity is 
extremely unlikely from blast loadirig, although 
it has been observed in long duration earthquake 
phenomena. Even if, for some reason, partial res- 
onance is achieved, the dampening in the system 


be~xcited will reduce the peaks considerably. 
_ 
In ,general, the combined effect of the 
two input motions, systematic and random, de- 
pends on 
their individual effects. It can be 
shown that the 
combined spectrum 
will be 
either equal to or less than the sum of the abso- 
lute values of the spectra corresponding to the 
individual inputs. It also appears reasonable that 
the combined spectre.rr: can be :xpected to be 
approximately equal to the square root of the 
sum of the squares of the individual spectra, 
point by point. In most practical cases of the 
type under consideration, the frequencies for 
which the spectrum values are important differ 
by a considerable amount, and the sums of the 
spectra or the square root of the' sums of the 
squares are nearly the same as the maximum in· 
dividual modal value. 


2-66 
Shock Spectra fo~ Field 
Ground Motion 
• 


• 
Using the concepts discussed above, it is 


po=re to derive shock spectra that can be used 
in design to assess the relative effects in a struc· 
ture (above· or below ground) or the effects on 
equipment within a structure. Such spectra can 
be described best by the use of a logarithmic 
plot 
that 
pennits 
values 
of displacement, 
pseudo-velocity, and acceleration to be read 
versus frequency in accordance with the rela- 


ti~ven by the equations in paragraph 2-65. 
_. In accordance with the limiting condi· 
tions discussed for Figure 2-89, the response 
spectrum can be represented by three regions, 
each region defmed by a straight line constitut- 
ing an envelope to the actual spectrum. An 


II 
example 
of such a spectrum for air-induced 
shock ior conditions of I00 psi surface pressure, 
5 
~ft yield. and an acoustic velocity of 2,500 
ft Isec is shown in Figure 2-90 for situLltions at 
the surface and at a depth' of) 00 f1. The left 
hand region. at the low frequency side, is influ- 
enced primarily by the maximum elastic trans- 
ient displacement. The intermediate part of the 
spectrum (horizon tal line) is a function of the, 
maximum velocity reached in the free field. The 
right hand side depends .on the maximum accel- 
eration. In brief. the approximate response spec- 
trum "envelope':' is described by three straight 
lines (actually there are discrepancies that may 
be as high as a factor of 2 in some areas. but the 
actual grolind motions are now known even thi, 
accurately. and' the simplification of the spec- 
trum in this way is permissible). 
• 
• 
• 
t 


• 
A line· D = constant, parallel to the dis- 
placement scales. drawn with a magnitude 
equal 
to the maximum ground displace- 
ment. 
• -A line 
V =constant, drawn with a mag- 
nitude of 1.5 times the maximum ground 
velocity. 


• 
A line A =constant, parallel to the accel- 
eration scales, drawn with a magnitude 
equal . to 
twice 
the 
maximum 
ground 
acceleration. 
_ 
Spectra may be sketched for horizontal 


an~tical motions at various depths and em- 
ployed as appropriate. Detailed procedures for 
constructing 
more 
accurate 
spectra 
may 
be 
found in "Effect of Inelastic Behaviot on the 
Response of Simple Systems. to Earthquake Mo- 
tions," A. S. Veletsos and N. M. Newmark, and 
"Response 
Spectra Approach to Behavior of 


Shock Isolation Systems," A. S. Veletsos and N. 
M. Newmark (see bibliography). In general, the 
air-induced and direct-transmitted ground mo- 
tions, as well as the out-running ground motions, 
if applicable, are computed. The spectra of each 
are plotted, and the design is made according to 
the envelop of the spectra. If it is known that 
the motions will arrive at significantly differen t 
times, the design or analysis is only made for the 
worst case. 
II The 
response 
spectrum. such as that 
shown in Figure 
2~90. may be used directly for 
elastic systems. If the frequency of equipment 
within a structure 'and its type and mode of sup- 
port are known, it is possible to ascertain certain 
results in terms of the response of the system. 
For 
example, 
for 
a 
single-degree-of-freedom 
system, the numper of gravitIes of ac~eleration 
corresponding to the natural frequencies may be 
considered as applied to the body to ascertain 
the dynamic str'esses. -For more complex sys- 
tems, the analysis is more difficult, although 
consideration of the most significant mode will 
lead to useful estimates of dynamic loadings in 


ma~cases. 


For. purposes of comparison and inter- 
est, 
ee shock spectra determined from reed 


- gages or ground motion records in atomic field 
tests are shown in Figure 2-91. The Operation 
HARDTACK spectra exhibit high accelerations; 
this is likely to be the result of the high water 
table that tends to emphasize the high frequency 
end of the spectrum. The maximum velocity is 
much higher than would be expected in a typical 
soil site. The Operation PLUMBBOB data ap- 
proximate the type of response that would be 
expected for more nonnal soil conditions. 


2-'1fi. 
Vs 


..., ....~.-:-. 


,t:i 
'ii 
>= 


~• 
..a:: 
'" 
.., 
.. .. 
~ - 
Q.- 
t 
W 
o~ 
·Vi 
:I 
Q.'tlc: 
0- 
o 
.:. 
-< 
:!!. 
I 


ti 
c: 
!. 
0 
en 
.~ 
o 
~....,o 


~~.;; 


~ t 
c: > 
8.- 
'" ... 
.. 
Cl. 
a:- 


.~ 


& 
>u 
Z 
I.IJ:lo 
I.IJ 
Ct 
II. 


;-. 
,... 
- 
~ 
~' 
,,,,A" 
.~- 


, 
, 
,". 
)0 
• 


t,,~!.~.·.). 
~ 
J . 
. ~.'.' .. ' ... ' .. ,~~ ..~ ... ~...~ 
I· 
;. 
• 
)0 
, 
• 
~ 
';..',' 
: • 


r 
~:-,.; 
y 
:. 
~~··.·..i- ./>'~<.>:v 
_... .'.r 
:.~ ~,.. 
,. 


• 
I 
( . 
' \. 
•• 
• 
•• 
: 
". 
". 
.. 
' .. 


l':.!.::,,"~,:, ...!.. .. 
~ 
.! •.- ...... 
:~~ 
~~ .• : ..I;'... 
.,-t.~~. ".::..'" 


h'" ,')' <Ii 
• 
, 
0{ 
'. t(. 
)0.") 
. 
~, 
'. '" 
>.. 
I( 
• 
oC 


. +'.:.>...:."\ 
.'~ 
., ......'..... 
. . ':./ 
:~:::'~<'-l .: 
' 
•. 
_ ............ '7' 
or· 
• '";~ ~ ;·f;::;.i.·.... 
.. 
.,.J •.,. ~.Y ..t 
..,1 
'. 
f· :..... .,. ,... .• •. A· or 
- 
.c( 
- 
u.·-"""f. "rJ" 
"'" 
. _.,...( •..~ 
~ 
"'" 
'.' _ 
.4 
I"" 
) .. ~.,. 
,.,.... 
.• 
.)if, r.; 
~ ,.. 
~ I.,A , 
. 


V .J... 
.! ..-' .!......... ."J. 
"" 
...1.1 ) ... j vi 
,.t.,.A. 
A 
."U ) ... .t 
.J. ".4' • 
). 
"'" 


• 
I 
;. 
• 
.. 
\ 
.. 
q. .,' 
I~' 
..)io 
.) 
.. 
\ 
.. , 
, 
} 
• 
. 


0(' ;'; ~ 
~. 
r 
.,:.- :~. ,.~~< . r.;~: 
11( .... <; "':-- . 
'.--: :'T.:". ~ .•'..~ .. ::.- 
~. '.:': 


~. /. 
. 
~ 
..... "'. ry • .)..! )0-. '" 
-.' 
•• 
'" 
..... 
. • -".'" ...... 
... ,. 
~ ... - 
.~••• 
), 
.. 
,,~'.... 
~ 
.. 
I 
.. 
• 
• 
". I.. 
!> 
• 
~ 
.' 
• "1" " .. 


... .!... oe.. 
..I • 
.!. 
...·:·~·.;,.."J1)..~.. j.;" 
)~.~.;. 
)yL 
i. ). 
~,.:.,,~..... 


I 
, 
.' 
• 
'....... 
• 
• 
• >. 
. . '..... 
.0 
~. ,0(" 
, 
,)r 
,)0 
" 
,,,...' 
I{..r... 
• 
r 
,.. 
f 
0(' 
, 
,. 
• fiT 
" 
• 
. 


• 
.. 
• 
,.' I 
; 
~. 
0·'''''''' 
.' 
• 
'. 
"'. ' 
.'. 
.J' '}. " 
~~o(I 
.. 
.,) ,'•• - .....'. -,.&- : •• 
~ 
,.. 
'" 
....'.... 
• 
.'. 
~. 
: .....". 
~... 
'. 


:'"~ 
.:~ 
'.~ 
".' 
"( 
'~:'.~ 
'\ 
",' 
"""t( 
." 
' .. 
"": 
".-"~ 
" 
..J.. 
~. ~ . 
.1... _ - 
...~•. ..,)... .Jw.~. } •• ". 
·.. ~·l 
.. .J•• ,.. . 
.J..... 
, 
.. 
, 
,.'." 
. 
" 
,).. >'.'""" 
." 
r , ("' 
~ ••'". 
" 


''b: '\. 
, 
.. "" 0( 
• 
I 
( 
,'I"..""" ( . , 
C 
.' 
• 
,.'." • 
• 
• 


,l-:':...'~ 
~...;....';..•..;....''''.:--.---·T--......p;...,'r·~---'...,''15""'--'--'r-----"l!'l""-'~·_''"-;;;_"-)-·-·_;·,-·--· _-:- 
..J" 


i 
•J 


, 
:! 


Iil;l 
.'j 


2-T8'S 
21(, 


,,., 
II 
I 
• 
(:l9S!,U!1 AJ.1::10131\ 
I I 
I 
• 
I 


Ee- 


iiju'c, 
>.... 
II 


2-'188. 
ZI7 


Problem 2-27. 
Calculation of the Ground Motions from 
a Nuclear Explosion at the Surface 


j. 


II The, nature of the free field effects of 
ground motion are described in paragraphs 2-60 
through 2-66. The attenuation 
fa~tor Q that is 
defined in paragraph 2-62 is used in the example 
presented below. As noted in paragraph 2-62. it 
generally is preferable to use the factor 1/(3 for 
attenuation; where (3 is described in Section II of .. 
Chapter 
11, 
and more specifically in Figure 
11-29. However, 1/(3 attenuation is not compati- 
ble with the equipment vulnerability curves in 
Figures 11-48 through II-50. Therefore, the fol- 
lowing example illustrates the uSe of the attenu- 


--ation factor described in paragraph 2-62; use of 
the 
1/(3 attenuation is illustrated in Problems 


IIJ-and 11,-4. II" 
. 


• 
Example 
.. 


Given: A I Mt weapon burston the surface 
of a medium that has a seismic velocity of 5,000 


.. 
. 


fl/sec. 
Find: ,The vertical and horizontal displace- , 
ment, 
velo.~ity,. and acceleration at a point 50 
feet below the surface and at a ground distance 
of 1,500 feet from the explosion. 
Solution: The equivalel)t distance for 1 kt 
. 
.. 


.=9 (1 ,000\°.5 (I ,000)(J)1/3 = 5.7 inches. 


.. 
100) 
5,000 


.. 
. 


It is reasonable to assume that impulse is pre- 
served in the stiff medium deScribed in this case. 


.. Therefore, the vertical elastic displacement at a 
depth of 50 feet will be the same as that at the 
surface, i.e., 5.7 inches. As described in para- 
graph 2-62, the permanent displacement at a 
depth of SO feet is 


d 
=Pao - 40 (1,000)2 (100 - y-) 
yp 
.. 
30 
c 
100 
p 
, 


=(1,000 .:. 40\ ( I,000)2 (100 - 50) 
30 
) 
5,000 
100 


= 0.6 inches. 


Thus, from paragraph 2-62 the total vertical dis- 
placement is 


dy = drJe + d yp = 5.7 + 0.6 = 6.3 inches. 


From Figure 2-17, the overpressure at the sur- 
face is about 1,000 psi. 
a.Air Induced Motions: The maximum 


elastic vertical displacement at the ground sur- 
face is given by the equation shown in paragraph 
2-62: 


C 
,., 9 (~IO ,.5 (1,000 )(WCMt))ll3 


Ie 
\100 j 
cp 


I 
i.. 


d 
(W(kt») 1/3 
= 
1,500 


(1,000)1/3 = 1SO feet. 


The peak vertical velocity given in paragraph 
2·62 is: 


where 


1 
Q = -....;;...---" 
y' 
+-H 


2-'t9Q. 
21B 


• 
and 


H =' 230 ( ~~ r.s 
(W(Mt)) 1/3. 


This gives 


H = (230) ,( I ~ggo r.s 
(1 )1/3 = 73 


1 


O! =--__ = 0.59 
I' 
50 
+'73 


v 
= 50 (1.000) ( 1,000) = 100 in./sec, 
5 
lOa 
5,000 


Vy = (0.59)(100)= 59 in./sec. 


The vertical acceleration is given by the equation 
in paragraph 2-62 as: 


ay =5C~~)(I~O) a 


= (5) (1,000) (100) (0.59) = 59 g. 
100 
50 


As described in paragraph 2-62, the horizontal 
values of the motions are taken to be 1/3, 2/3, 
and I times the vertical displacement, velocity 
and acceleration, respectively, ioe., 


d 
= 1/3 d 
= (1.) (6.3)= 2.1 
inches, 
x 
y, 
3 
' 
. 


Vx = ; vy = (;) (59) = 39 in./sec, 


b. 
Direct Transmitted Motions: The ex- 
pressions given in paragraph 2-64 are applicable 
for direct transmitted motions. The radial mo- 
tions are essentially horizontal, for the distance 
and depth of this.example. These are (note the 
distance is less than 2,500 ft/Mt 1/3 ): 


( 
)1.5' 
d 
I = 4 (W(Mt») 5/6 
I,O~O 
. 


;., (4) (1)5/6 
(.1,000)1.5 = 22 inches. 
. 
1,500 
0 


VI 
=.12 (W(Mt») 5/6' (l,O~oy.s 
C.~~O) 


= 12 (1)5/6 (1,000)2.5 (5,000 ), 


1,500 
'1,000, 
= 22 in./sec, 
, 


ar =' 0.4 (~(Mt)) 516 C,o~oy.5 C 
,~oor 


= (0.4) (1)5/6 (1.000)3.5 (5.000)2 
1,500 
1,000 


= 2.4 g. 


In this case, the motions normal to the radial 
motions may be taken to be vertical. These are 
1/3, 2/3, and 1 times the corresponding radial 
values of displacement, velocity and accele,ra- 
tion, i.e., . 


dy =+dr 
==(n (2.2) = 0.7 inches, 


vy = ; 
vr = (; )(22) = 15 ino/sec, 


2-...... 
Z/' 


3,' 


• Answer:· The peak values of· both hori- 


zontal and vertical motions, for air-induced and 
direct-transmitted shocks are: 


d(in.) 
v(in./sec) 
~ 


Vertical· 


Air-induced 
6.3 
59 
59 
pirect-transmitted· 
0.7 
15 
2.4 


Horizontal 


Air-induce'd 
2.1 
39 
59 
Direct-transmitted 
2.2 
22 
2.4 


( 
. In this particular example the air-induced mo- 
tions; with the.exceptivn of the horizontal dis- 
placement, are the largest values for both the 
vertical and 
i horizontal directions. In the hori- 


--'zontal direction, the direct-transmitted value of 
displacement IS nearly equal to that of the air-· 
induced 
valu~. In many cases, especially at short 
range and for stiff media; it will be found that 
the direct-transmitted motions will be the larg- 
est, and will control the analysis. From para- 
graph 2-66, the bounds for the shock response 


spectra are obtained .by multiplying the cOntrol- 
ling 
displacement, 
velocity. and 
acceleration 
values by 1, LS, an:! 2.0, respectively. 


VertiCal Response Spectrum Bounds: 
D = I.Od = 6.3 in. 
Y = 1.5v = 89 in./sec 


A = 2.<kl = 118g 


Horizontal Response Spectrum Bounds: 
D = l.Od = 2.2 in. 
V = 1.5v = 59 in./sec 


A = 2.<kl = 118g 


These spectra bounds are shown in Figure II-51. 
They show the response spectra for a single- 


degree~f-freedom 
system. These 
bounds are 
used in the analysis of the vulnerability of equip- 
ment mounted in an underground structure in 
Problem 11-8. 
• 
Related Material: See paragraphs 2-24, 


2.~hrough 2-66, and· Sections II and III of 
Chapter 11. 


- . 
SECTION IV 


_ 
UNDERWATER EXPLOSION 


PHENOMENA '. 


• 
An underw;iter nuclearexplosion is char- 
acterized initially by the formation and .propaga- 
tion of an underwa'ter shock wave. The dissipa- 
tion of heat at the shock front causes the water 
to vaporize in the vicinity of the explosion and 
results in the fonnation of' a large bubble of 
steam. The interaction of the bubble and the 
shockwave with the air-water interface produces 
various 
surface 
phenomena, 
including 
spray 
domes. columns, plumes, base surges, radioactive 
pools and surface waves. If the explosion is on 
o~r the bottom, a crater is fonned. 
_ 
An underwater nuclear explosion also 
can result in thennal, nuclear, and electromag- 
netic radiation above the surface, but these phe- 
nomena are insignificant, except in the case of a 
shallow burst. These phenomena, when caused 
by a shallow burst, should,approximate, those 
from a shallow underground explosion. 


2-67 
Underwater Bubble III 
II An 
underwater 
nuclear explosion 
~­ 


leases large amounts of thennal and nuclear radi- 
ation, essentially all of which is absorbed by the 
surrounding water within several feet of the ex- 
plosion (some radiation in the visible spectrum 
can be radiated to greater distances, depending 
,olJ.lii transparency of the water). 
, • 
During the early stages of the explosion, 
the warhead materials attain a very high temper- 
ature (tens of millions of degrees) and a very 
high pressure (on the, order of millions of atmo- 
spheres). Energy acquired 'by these materials is 
transferred' to the layer of water nearest the 
bomb, 
which is heated and compressed and 
which, then, heats and compresses the next out- 
. ward layer. A compression wave (or hydrody- 
,narnic shock front) is formed by this mechan- 
ism, and this wave moves outward from the 


explosion at a speed greater than that of sound. 
This shock front expands faster than the ma- 
terial that.it engulfs, which also moves outward, 
b"a slower rate. 
• 
As the shock front moves away from the 
pomt of explosion, energy is dissipated as heat, 
which' raises the temperature of the water en- 
gulfed by the, shock front. The largest tempera- 
ture increase occurs near the center of the explo- 
sion,where the water is not only vaporized but 


dissociat~das well. At greater distances, the 
water is vaporized and turned to steam at the 
shock front., At still greater distances, the water 
is not heated by the shock front sufficientiy to 
boiL'The bubble that is formed does not have a 
clearly defined "surface:' Figure 2-92 shows this 
prliession in a highly stylized form. 
Thus, shortly after an underwater burst 
an expanding bubble ,is formed, which is com- 
posed largely of vaporized water with radio- 
active bomb debris at its center, surrounded by 
heated 
water. 
Large 
quantities 
of 
neutron- 
activated elements in the sea water are also pres- 
ent in the steam and layer of surrounding water; 
however, the total induced radioactivity is rela- 
tively small compared to the fission product 
radioactivity. 
Continued 
expansion 
of 
this' 
bubble reduces the pressure within it. As the 
bubble pressure falls below the vapor pressure of 
the heated water,. additional water flashes to 
steam at the interface of the bubble and the 
vii If~he explosion is deep underwater, the 
bu 
e continues to expand at a decreasing rate 
until a maximum size is reached. Ifit is not too 
near the surfaCe or the bottom, the bubble re- 
mains nearly spherical to this pomt. !-oS a result 
of the inertia of the water surrounding the ex- 
panding bubble, the 
bubble actually overex- 
pands; i.e., when it attains.its maximum size, its 
pressure is ,below the ambient water pressure. 
The high pressure around the bubble then causes 
it to contract with an increase in its pressure and 
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Figure 2-92.II Representation of Bubble DevelopmentII 
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condensation of some of its steam. Since the 
hydrostatic pressure is larger at the bottom of 
the bubble than at the top. the bubble does not· 
remain spherical during the contracting phase. 
The bottom moves upward faster than the top 
(whic;h may even remain stationary) and con- 
tacts the top to form a toroidal, or doughnut- 
shaped. 
bubble 
as 
viewed 
from 
above. This 
causes t_urbulence and mixing of the bubble con- 
tents with 'the surrounding water.' 
.. 
_ 
The inertia of the water set in motion 'by 


th~bble contraction causes it to overcontract, 
and 
its internal 
pressure 
becomes very much 
higher than the ambient water pressure. A 
sec~ 


ond compression wave (the first bubble pulse) 


, commences when the' bubble reaches its mini- 
mum volume. This wave has a lower peak pres- 
sure. 
bu t a longer dura tion. than 
the 
initial 
shock. A second cycle of bubble expansion and 
con traction begins. 


• 
When the bubble pulse is emitted, the 


bu= surface becomes exceedingly unstable; a 
multitude 
of water spray jets shoot inward, 
thereby condensing much of the steam 'and 
de~ 


grading the bubble energy. Figure 2-9.3 shows 
the relative partitioning of the energy of the 


ex~osion with time qualitatively. 
The bubble is relatively stationary dur- 
ing 
e initial expansion cycle, but upon con- 
tracting.it begins to migrate upward uncler the 
action of buoyant forces. The rate of upward' 
migration is greatest when. it is contracted and 
nel zero when it is expanded. 
,,' 


If the bubble is far enough below the 
surface, it continues to pulsate and rise, though 
after three complete cycles enough steam will 
have condensed to make additional pulsations 
unlikely. During the pulsation and upward mlr 
tion of the bubble, the water surrounding the 
bubble acquires considerable upward momen- 
tum and eventually breaks through the surface 
with 
a 
high 
velocity 
(about' 200 
mph 
for 
Wigwam, a deep underwater shot) thereby creat- 


in~large plume. 
j 


If the explosion is at a shallow depth. 
the 
ubble may break through the surface dur- 
ing one of the early pulsations, or even during 
the fIrst pulsation cycle. If such a breakthrough 
occurs when the bubble pressure is higher than 
am bient (as with a very shallow explosion), a 
'phenomenon known as a "blowout", occurs. If a 
breakthrough occurs when 
bubble pressure. is 
below ambient, "blowin" occurs. The surface 
effects differ markedly for the, two phenomena 
(see Figures 2-100 and 2-101 and accompanying 
text). A migrating bubble will not produce a' 
bubble, pulse if the bubble minimum occ,urs at or " 
a~ the original water surface. 
, 


• 
A strongly migrating bubble might con- 
tinue to expand and not pass through a maxi- 
m urn 
size, 
somewhat 
like 
a 
balloon 
rising 
through the atmosphere. This is the balloon ef- 
fect. A bubble migrating in this manner will not 
emit' any 
bubble pulses. The conditions that 
cause a, nuclear 'bubble' to balloon rather than to 
pulsate are still subject to speculation: however, 
the absence of an expected bubble pulse might 
beKxlained by the balloon effect. 
'The maximum 
bubble 
radius and the 


fIrst 
ubble period (time between burst and first 
minimum) are related to yield and depth of 
burst, as illustrated by the following twoequa- 
tions, which pertain to submerged bursts, i.e., 
the depth of burst, db is greater than A 
: 
m a x, 


Wl/3 , 
A 
= J'-- 
max 
Z1/3 ' 


wl/ 3 
T. = K'-- 
Z5/6 ' 


where 


A~ax = Maximum bubble radius (ft), 
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Figure 2-93. II 
Energy Distribution Following an Underwater 


Nuclear Burst til 
' 


- 
J' = Empirical constant, 1500, 


K' = Empirical constant, 515, 


T) 
= First bubble period (sec), 


J-l,',= Yield (kt), 


Z = Hydrostatic head 
(db + 33) (ft). 


The form of the equations is the same as those 
given by,Cole for high explosives. The constants, 
J' 
and K', 
are 
slightly' different because an 
underwater nuclear burst forms a steam bubble 
rather than gas bubble (as do high explosives) 
and because of other differences, mainly the dis- 
tribution of total energy between bubble energy, 
shock 
wave 
energy, 
and 
residual 
radioactive 
enlr', 
' 
It is likely that the constants J' and K' 
are actually somewhat dependent on depth. This 
can be understood qualitatively as follows: J' 
and K' implicitly include the fraction of explo- 
sion energy includea in the bubble. As the depth 
of burst is increased, a smaller fraction of explo- 
sion energy is converted to steam by the expand- 
ing shock wave (and more energy goes into heat- 
ing water) because of the greater hydrostatic 
pressure that must be overcome to boil water. 
• 
The second and third bubble maximum 


ra~he periods between minima, and the mi- 
gration between periods, depend on the fraction 
of bubb'le energy that remains following suc- 
ceeding bubble minima. The energy -loss that 
occurs is due primarily to ste,am condensation at 
bubble minima as previously described, but the 
amount of energy loss is irifluenced by the depth 
at which the minm,um occurs. An analytical rep- 
resentation of the energy as a function of bubble 
migration, ~which is best solved by use of a high- 
speed computer, has been developed," however, 
satisfactory solutions to bubble problems can be 
obtained by use of curves derived from the rep- 


resentation, as illustrated in Problem 2-28. 
'I' 
WATER SHOCK WAVE AND OTHER 
P ESSU,RE PULSES II 
' 
II The initial shock wave from an under- 
water explosion propagates radially from the 
source. This shock wave is characterized by an 
abrupt, virtually instantaneous, increase in pres- 
sure fonowed 
by a decrease that is approxi- 
mately 
exponential. 
Near 
the explosion, the 
peak shock wave pressure is extremely high, but 
energy losses to the water cause the pressure to 
decrease somewhat faster than inversely with 
radius. Sim;!arly, ear!:' shock wave propagatic,'i1 
velocities are high, but by the time the peak 
pressure falls to about 3,000 psi, thepropaga- 
tion velocity becomes nearly equal to the acou- 
stic speed in water (about 5,000 ft/sec), and the 
enlrlosses become small. 
, , 
If an explosion takes place far from 
eit er the surface or the bottom, the shock wave 
remains spherical as long as it propagates into 
water having constant acoustic velocity. This ve- 
locity depends on temperature, hydrostatic pres- 
sure, and salinity, however, and in regions of the 
ocean in which significant sound velocity gra- 
dients exist, the shock wave can be bent or re- 
fracted, Refraction can either increase or de- 
crease the shock wave pressures locally, depend- 
ing on ambient water conditions. No generaJ pre- 
diction techniques can be given, but the subject 
is.ussed further in paragraph 2-71. 
' 
When a shock wave moving in water en- 
counters another medium, it may be reflected as 
a tensile or rarefaction wave, as in the case of a 
water-air 
boundary; a compression or shock 
wave, as in the case of a water-bottom interface 
if the angle of incidence is not too far from the 
normal; or a distorted pulse having in general 


;a;:;:-;-xpreuiOn iI cIeIc:ribed in ''The Parameters of Under- 
wllr;i;c~ Explolioa Bubbles" <_ bibliolraPhy). 
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II 
both a positive and a negative phase; as in the 
case of more glancing incidence at the bottom. 
These are discussed separately in the foUowing , 
paragraphs. 


2·68 
Surface Reflection II 
. 
It The rarefaction wave, generated by the 
reflection of the priinary shock wave from the 
surface, propagates downward and relieves the 
pressure behind the primary shock wave. If the 
shock wave ;is treated as a weak (acoustic) wave, 
this .interaction instantaneously decreases the 
pressure in the primary shock wave to a value 
that might be well below the vapor pressure of 
the water, as shown by the broken line in Figure 
2·94, Point A. Cavitation occurs in seawater 


--when its pressure decreases to a value somewhat 
below its vapor pressure. The pressure of the 
primary shock wave is, therefore, reduced to a 
value which, when compared with the peak pres- 
sure, is usually so close to ambient water pres- 
sure that the shock wave pulse appears to have . 
been 
truncated, 
i.e., 
reduced 
to 
ambient 
pressure. 
I 


FaT a strong primary shock wave, the 
re ected rarefaction wave propagates into water 
that has alieady been set in motion by the shock 
wave, Therefore, the rarefaction wave arrives 
.:....;:;;;; 
~~~::-. pr~dicted from the acoustic approxi- 
mation, and the pressure cutoff is not instan- 
taneous. This effect typically gives a pulse shape 
shown by the solid line for Point A of Figure 
2-94. The shallower the point at which pressure 
measurements are made, the sooner the primary 
shock pulse is "cut ofr' and, hence,the shorter 
its duration (see Figure .2-94, Point B). At suffi· 
ciently shallow locations, the rarefaction wave 
interacts with the shock front and reduces the 
peak pressure (see Figure 2·94, Points C and D). 
The region in which peak pressure is reduced is 
known as the uanomalous region." 
. 
I 


TIle effects of surface reflection decrease 
rap 
y with increased depth of either the explo- 


sian or the point of measurement. Conversely, as 
the depth of burst is decreased (or the yield in- 
creased for a given depth of burst), the effects 
increase. The size of the anomalous region in- 
creases with decreased depth of burst u~!!1' for a 
surface burst, the anomalous region ::9fuU, all 
points beneath the water surface except those 
close to the explosion and directly under it.· 
tI A limited amount of data are available 
concerning the reflection of shock waves from 
an ice layer. To date all tests have been with 
relatively small explosive charges (a majority of 
the test shots have been in the)· to 4D-lb range, 
with one shot of 630 Ib).The situation is com- 
plicated by reflection and refractions E.: both the 
water-ice boundary and the ice-air boundary. 


2-69 
Bottom Reflection • 


__ Under certain circumstances, the shock 
wave from the bottom can be more damaging to 
surface ships and shallow submarines than the 
primary shock. Although the peak pressure of 
the reflected wave is usually smaller, the wave 
arrives at a steeper angle, and therefore may in- 
duce more damaging shock motion to a target at 
or near the surface. 


2·70 
Secondary Shocta.lnd 
Pressure WIVeS • 


__ 
An underwater burst can cause compres- 
sion waves in addition to those described above. 
but these effects are usually negligible. These 
waves include retransmitted pulses, cavitation 
pulses, bubble pulses, and others that are de- 
scribed briefly below. 
• 
Retransmitled Pulses. 
Upon 
reflecting 
from the surface and bottom, some pri- 
mary shock energy is transmitted into the 
air or bottom material,and some of this 
energy can be retransmitted back into the 
water. 
• 
OlvitaNo" Effects. AI noted earlier, the 
reflection of the shock wave from the sur· 
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II face can lead to bulk cavitation below the 


surface, During the period of cavitation, 
the cavitated region can absorb funher 
shock waves impinging on it. In 'the forma- 
tion ofl the cavitation region, the rupture 
of the ;water by the shock wave causes a 
liquid 
surface 
layer 
(sometimes 
called 


"spall") to project upward in almost a bal·' 
listic trajectory. The impact of this spall 
when it returns to the surface can lead to 
further 
secondary 
shock 
waves. 
These 
shock 'waves account for the occasiohal 
larger 'damage 
to 
ships at intermediate 
I 
ranges 'compared with the damage at short- 
er or longer ranges. 
• 
Bubble Pulses. 
If an explosion is deep 
enough for one or more bubble pulsations 
to occ~r, compression waves are generated 
at the time of ~ach bubble minimum., 


• 
Other Pulses. All pressure waves are sub- 
ject to multiple reflection from the surface 
and bottom. Also, in areasof irregular bot- 
tom topography, more than one reflected 
pulse Can be generated: 


_ 
Figure 2-95 shows typical shock wave 
and pressure pulse patterns for various burst and' 
measurement conditions. 


2.~ Refraction of Shock Waves _ 


• 
The shock front from an explosion far 
from a boundary (surface, bottom, etc.) remains 
spherical provided that it expands into water 
having constant acoustic velocity (isovelocity). 
However, the acoustic velocity in water depends 
on the temperature, salinity, and hydrostatic 
pressure, and these properties are not necessarily 
constant throughout a large body of water. As a 
result of variations in these properties, generally 
in horizontal layers, a region of water can have a 
characteristic acoustic velocity proflle (sound ve- 
locity vs depth). Under such conditions, the 
shock wave is bent (refracted) because one part 
of the wave moves faster than another. 


_ 
Refraction of the shock wave can result 
in convergence and reinforcement of the shock 
wave. \ This 
reinforcement 
commonly 
occurs 


along a surface called a "caustic." This effect 
may be illustrated by the use of a ray diagram· 
(Figure 2·96). The acoustic velocity prome that 
corresponds to the ray diagram is also shown in , 
the fIgUre. The formation of several caustics for 
one burst is observed frequently. Regions of 
relatively low pressure also occur, such as the 
"shadow zone" above the caustic in the fIgUre. 


• 
The caustic can reinforce the peak pres- 


sur~y a factor of five or more over the isove- 
tocity value, while the pressure observed deep in 
the shadow zone can be essentially negligible. 
Shock wave impulse amplification factors do not 
go through a maximum at the caustic position. 
They remain much closer to unity than do peak 
pressure and energy factors. 


__ As a result of ocean currents and under- 
water swells, the acoustic velocity proflle in the 
ocean is rarely stable: Thus, the location of the 
caustics can shift rather quickly and unpredict- 


ab~similar.to optic mirages). 
• 
Common types of sound velocity pre; 
flles can give rise to caustics at various distances 
greater than about one mile. In particular, in 
about half of the areas of the oceans, a caustic 
occurs at what is known as the convergence 
zone, typically 30 (10 to 40) miles from an 
underwater burst. Therefore, it is conceivable 
that a submarine might suffer damage from its 
own multi-kiloton weapon detonated at a dis- 
tance of about 30 miles if it happened to lie on 
the caustic of its shot. 
II In Figure 2·96, sound speed is consider- 
ed to vary with depth only. Actually this would 
. rarely happen. The velocity can vary consider- 
ably across the ocean within the nnge of inter· 
est. Generally, the region of focusing is sharper 
~ .-.ted by • CIIIrW abat II atny. perpeD4k:uIIr 
tollt~~ ahock froGt. 
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Figure 2-95. II Typical Shock Wave Patterns 'Along Line A-A 
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• 
and the maximum value of the amplification is 
greater for smaller yields. 


2~72 
Air Blast II' 
. 
. 


. 
• 
The 
air blast wave has one to 
three 


pu~of various amplitudes, depending on the 
.burst depth and the location of the observation 
point 
with 
reference 
to 
surface zero. These 
pulses are generated by two .underwater explo- 
sion mechanisms: One is the underwater shock 
wave 
that 
transmits a portion of its energy 
across the water-air interface. The other is the 
underwater explosion bubble, which manifests 
itself in two ways: for burst depths shallower 
than about 
351~'ll3 feet. the bubble vents. caus- 
ing an ai, pressure jJulse: for Durst dept.is deeper 
than 3511'113 feet. a spray dome is pushed up 9Y 
tJxpanding bubble, causing a bow wave. 
For burst depths shallower than about 
50 
3 feet the transmitted pulse arid the bub- 
ble generated pulse are formed almost simultan- 
eously. 
resulting 
in 
a single intense pressure 


. pulse over surface zero. For burst depths be- 
.tween 5011'1/3 feet and about 15QWl13 feet, up 
to 
three 
pulses can -be seen. The transmitted 
shock in air bifurcates near surface zero to fonn 
two air pulses.' and the bow shock causes the 
third air pulse. For bUTst depths. deeper than 
about 
15011'1 /3 
feet, the spray dome rises so 
slowly that no bow shock is formed and only the 
air pulse from the transmitted shock is important. 


• 
This discussion of the generation of air 


bla~rom underwater bursts is a simplified dis- 
cussion of a complex subject. Much ,analytical, 
theoretical, 'and· experimental work has been 
done to understand the physics of air blast gen- 
eration. More complete discussions of the sub- 
ject are given in Malme, Carbonell and Dyer 
(1966); Peckham and Pittman (1968); Pittman 
(1968); and Rudlin and Silva (1960) (see biblio- 
graphy). Only a limited number of air blast 
measurements are available from nuclear weapon 
tests. The data are not sufficient to fonn a basis 


for predicting air blast for a spectrum of burst 
depths and pressure levels. The available pte- 
,dictive measures are based mainly on high explo- 
"sive tests as described in Problem 2-31. 
II 
SURFACE EFFECTS OTHER 
THAN WAVES • 
' 
II Thenrst surface effect of an underwater 
burst is caused. by the intersection of the pri- 
mary shock wave and the surface. Viewed from 
above, the effect frequently appears to be a 
rapidly expanding ring of darkened water, often 
called the "slick" (Figure 2-97). A white circular 
patch (the "crack") follows closely behind the 
darkened region. The crack probably is caused 
by underwater cavitation produced by the rc- 
fleeted rarefaction wave. Shortly after appear- 
ance of the crack. the water above the explosion 
rises vertically·. and forms a white mound of 
spray, called the "spray dome" (Figure 2-98). 
This dome is caused by the velocity imparted to 
the water near the surface by the ret1ection of 
the, shock wave and the subsequent breakup of 
the surface layer into drops of spray. Additional 
slick, crack, and spray dome phenomena may 
result if the. bottom-reflected shock waves and 
bubble-pulse compression waves reach the- sur- 
- fawith sufficient intensity (Figure 2-99). 
_ 
The spray dome from a shallow burst 
changes rapidly to a shallow column formed by 
the upward and outward motion' of the water 
surrounding the explosion bubble.. If blowout 
occurs, the upper part of the column is likely to 
resemble a crown, which contains explosion pro- 
ducts blown out of the column (Figure 2-1001. 
If blowin occurs, the crown is likely to be absent 
(Figure 2-101 a). In its later stages. the colur.1n 
may break up into column jets (relatively broad 
spouts of water that disintegrate into sprJy 
~\S 


thMtravel through the air) (Figure 2·10 Ib). 
If the burst is sufficiently deep to pre- 
vent 
lowout, but shallow enough for the bub- 
ble to continue to oscillate as it approaches the 
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"Figure 2':100. _ 
Very .Shallow Burst Showing Crown of Explosion 
on Column (oblique aerial view 0.7. sec after burst) III 
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Figure' 2-101. II Column from Shallow Underwater Bunt (no blowout) II 
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Table 2-8. 
_ 
Classification of Bunt Depth 
Categories II 


"'&hese values are derived from H.E. tests. They have 
nKen confirmed by nuclear tests; therefore, they are 
to be interpreted as "containment possible" and "possIble 
limiting depth for spray dome formation." 


hazards. 
• 
Very 
shallow bursts. 
The bubble vents 
early during the first expansion cycle, i.e., 
when the bubble pressure is greater than 
ambient, and fission products are' blown 
out at that time. 
• 
Shallow 
bursts. 
The bubble 
vents 
latt:: 


during the first cycle,after the bubble 
pressure has dropped below ambient, and 
fission products are blown in. 
• 
Deep bursts. The depth of burst is equal to 
or greater than the radius of the fully ex- 
panded bubble, but not as deep as the very 
deep burst described below. 
• 
Very deep burst. The explosion is at suffi- 


o cient depth that the bubble breaks up, be- 
comes a vortex ring, or loses its identity' 


-£efore reaching the surface. 
~ The surface radioactivity from a very 
deep burst is operationally insignificant within 
about an hour and may be faint enough to es- 
cape detection. 


II 
surface. plumes of upwelled water and steam can 
occur (Figure 2-102). 


_ 
Upon 
subsidence 
of the column and 
plumes 
from 
an 
underwater 
explosion, 
a 
"doughnut-shaped 
ring" 
of 
mist 
is 
formed, 
which is called the base surge and which is high- 
ly radioactive (Figure 2-103). This ring, or series 
of rings, expands radially in the absence of wind, 
but in the ,presence of wind, it elongates and 
drifts with the wind until it dissipates. A train of 
surface waves also expands radially from the ex- 
plosion. 


_ 
After dissipation of the base surge, the 
water 
surface 
around 
the 
explosion 
appears 


white (Figure 2-104). This "foam-patch" results 
from the upward motion of water and. gaseous 
explosion products in the vicinity of the bubble, 
their spreading over the surface of the patch, 
and their downward motion at the edge of the 
patch. During its later stages, the foam-patch 
appears as a ring of foam and debris that is left 
floating where the water circulates downward. 


_ 
As previously noted, the surface phe- 
nomena described above vary with the weapon 
yield arid the depth of burst. Although clear-cut 
distinctions cannot be made in all cases, six ex- 
plosion categories have been established to aid in 
establishing safe w~ap6n delivery criteria. These 
categories are defined in Table 2-8. The classifi- 
cations have been tested at yields as large as 30 
kt, but they probably are valid for yields up to 
about 100 kt. 
• 
The limits of the relations given in Table 
2-8 are plotted in Figure 2-1 OS, and the phe- 
nomena that establish the categories are describ- 
ed below. 
• 
Near-surface bursts. 
The layer of water 
above the burst is vaporized by the explo- 
sion. The surface phenomena for this type 
of burst and the associated hazards are 
unknown. The radiological hazard of the 
base surge is considered to be unimportant 
compared 
with 
air 
blast 
and 
fallout 


Classification , 
00 


Near surface 


Very shallow 


Shallow 


Deep 


Very deep 


Contained'" 


No spray dome'" 


DeflIlition 


113 
o < db < 
21 W 
21WI/3 < d < 
75WI/3 
b 
75WI/3 < d < 24OWI/4 
b 
240WI/4 < d < 7ooWI/4 
b 
7ooWI/4 < db 
1550Wl/4 < db 
25OOWI/3 < db 
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Fi~re 2-102. 
_ 
Plume from 'Upwelling of Water from Very Deep 
Burst (24 sec after· burst) 
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Figure 2-105. 
_ 
Classification of Underwater, Nuclear Bursts • 
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• 
No radioactive debris would be observed 
on the surface from contained bursts. The prob- 
ability of containment of the radioactive debris 
is a function of the oceanic density prome. The 
relation in Figure 2-105 probably is valid in cer- 
tain tropical, waters, but bursts at these same 
depths in ,northerly latitudes probably would 
nu contained. 
_ 
No ronace effects due to bubble pulsa- 
tion or migration occur as a result of bursts that 
do not produce a spray dome. 
' 
• 
Although water depths may also influ- 
en ce 
surface 
phenomena, 
and 
consequently 
burst classification, the mathematical defmitions 


in~le 2-8 do not include this influence. 


_ 
Thevarlous surface phenomena can be 
__divided into: two categories aCCQrdingto their 
origin, i.e., those produced by shock waves and 
those 
produce~ by the water motion accom- 
panying bubble pulsation and rise. In the former 
category, the slick, crack; and spray dome are of 
operational importance. In the latter eategory, 
columns, plumes, base surge, and the long-lasting 
foam patch can be significant. More detaileli de- 
scriptions of the surface phenomena are given in 
the following paragraphs. 


2-73 
Spray Dome 
_ 
'!! The spray dome (or spall), is formed by 
the 
acceleration 
and breakup of the surface 
water layer. The initial surface velocity directly 
above the burst is the sum of the water pl£rticle 
velocities imparted by the incident and reflected 
waves; elsewhere, the initial, velocity is approxi- 
mately the sum of the vertical components of 
these velocities (!'.ome initial surface velocities 
have been observed to be higher than so calcu- 
lated). The, height of the spray dome is limited 
by deceleration of gravity and air drag. 
_ 
If the depth of burst is sufficiently deep 


fo~e or more bubble oscillations to occur be- 
fo... venting, 
there may be secondary spray I 
dome-s from reflection of compression waves 


emitted by the later bubble expansions. These 
can be" several times higher than the primary 
dome from very deep bursts. 


• 
The calculations of spray dome extent, 
ve leal water velocity, and spray height are de- 
scribed in Problem 2-32. Plumes appear over the 
central portion of the spray dome from shallow 
bursts; therefore, a calculation of spray dome 
height may have little value. 


2·74 
Plumes, Column, 
_ 
' 
Cluliflower Cloud ~ 
II If the burst is deep enough to prevent' 
blowout (shallow to very deep bursts), expan- 
sion, pulsation, and rise of the bubble produce a 
violent upheaval and expansion of water above 
the surface to form plumes. In those deep explo- 
sions where the bubble collapses just beneath 
the surface, the rising bottom of the bubble jets 
through the layer of water above the bubble and 
probably will form a plume. In deep and very 
deep 
explosions, 
the 
upward motion of the 
water about 
the rising bubble also forms a 
PII' The water displaced by the expanding 
bu 
e from very shallow bursts forms a dome- 
shaped shell of water over the bubble. Later, this 
Shell ruptures and blows out the bubble constit- 
uents. For a very shallow nuclear explosion, 
blowout is accompahied by the- formation of a 
"cauliflower cloud" at the top of the plume. 
II 
Since there are insufficient data for a 


comprehensive prediction technique, predictions 
of plumes, columns, and clouds are based on 
direct scaling from four nuclear bursts. Pre- 
dictiverelationships and' examples of their use 
are given in Problem 2-33. 


2·75 
Bese Surge • 


• 
As the plume from a deep burst or the 
co urnn from a shallow burst falls back to the 
surface, it generates I 
cloud of fine spray, or 
lerosol, c:alled "base surge," .which usually is a 


o 


) 


II 
radiation hazard. The base surge expands and 
circulates as a toroidal cloud under the influence 
of the collapsing .vater plume. The expansion is 
rapid atftrst (exceeding 100 ft!sec); however, it 
decelerates upon mixing with the ambient air. 
During this period, part of the cloud evaporates 
and cools the air.' As its kinetic energy dissipates 
and the cloud reaches its maximum diameter, it 
tends to thin out, and it may even become invisi- 
ble ,as a consequence of diffusion and continued 
evaporation. 
In 
a sense, the base surge still 
exists, because some ocean salts, and detonation 
products remain airborne,. ' 
II Two types of base surge have been en- 
countered during weapon tests: 
• 
"Doughnut-shaped" clouds from very shal- 
low and shallow detonations 
• 
"Disc-shaped" clouds or, more properly, a 
series 
of . concentric 
"doughnut-shaped" 
clouds 
from deep 
aOnd very deep deto- 
_ 
ations. 
' 


Although 
visible 
configurations 
are 
actually quite irregular, the base surge configura- 
tion for very shallow and shallow detonations 
can be represented as a hollow cylinder'or a hol- 
low truncated cone Jor radiation exposure esti- 
mates. The inner diameter of the cylinder is ap- 
prOXimately twcHhirds the outer diameter. For 
very deep and deep bursts, a cylinder or trun- 
cated cone are adequate model conftgurations. 
"The height of the visible b~ surge clouds at 
various, weapon tests have varied between about 
1,000 and 2,000 ft, probably as a result of tur- 
bulence and meteorological conditions. 


_ 
A base surge is transported by the wind 
dunng its formation and growth. As a general 
rule, the transport should depend upon the wind 
speeds at altitudes between the surface and the 
maximum surge height. The downwind growth 
'curve may be predicted by adding the highest, 


wind speed reported between these levels to the 
average 
growth 
curves. 
The 
upwind 
growth 


should lie between a predicted curve based on 
subtracting half the surface wind speed (a rule- 
, of-thumb to 'account for friction with the sur- 
face) and a curve obtained by subtracting the 
wind speed used for the downwind curve. At 
early times the former would be more accurate 
while at' later times the latter would be more 
correct. Base surge prediction relationships and 
an example of their use are given in Problem 
,2-34. 
' 


2-76" 
Foam Patch and Ring III 


• 
After subsidence of the plume, column, 


ann'ase surge, a loughly circular area of dis- , 
turbed water is still visible. This area is primarily 
deftned by its white color and indications of tur- 
bulence at the perimeter. In the case of deep 
explosions, the patch can result from an upsurge 
of water caused by bubble migration. This water 
spreads out and then sinks, the overall motion 
being sOmewhat 
toroidal. Foam and 
flotsam 
tend to collect at the edge of the patch. 
, _ 
The whiteneSs from deep, shots is be- 


lie~ to be caused by foam generated by the 
agitation of the water, particularly as the plumes 
fall back. Similar foam patches occur from shal- 
low detonations, but their whiteness is caused 
partly by suspended bottom material. The foam 
patch contains appreciable quantities of weapon 
debris and presents a potential radiation hazard 
t0l:nnel in small boats or ships. 
, 


The foam patches developing from test 
exp osion have been on the order of 1.5 to 2.5 
miles in diameter about 15 min after detonation, 
, but they continued to expand slowly and irregu- 
larly with water currents and diffusion into the 
water, until the last visible manifestation was a 
ring of foam and debris. Their radioactivity 
could be traced for many days, although they 
had been,inVisible for quite some ,time. 
' 


~ 
WATER SURFACE WAVES 
_ 


2·n 'Generation and.oagation of W.18r 
Surface Waves 
' 
.. A submerged nuclear explosion, like any 
oM:"physical event that produces a localized 
disturbance on the water surface, generates a 
group of surface waves that expand radially. The 


cha.-act~:i:;ti=s of explosively generated waves 
depend upon the yield, the depth of burst, and 
the range to the point of 
observatio~. If the 
water is sufficiently shaHow, the bottom con- 
tour also can affect the generation and propaga- 
tion of the waves. 
II The 
following physical description of 
the wave generation' process and the succeeding 


__prediction technique are based on: 
• 
Laboratory experiments with small' H.E. 
charges 
• 
Wave data for the nine nuclear shots that 
have been fired on or in the lagoons at 
Eniwetok and Bikini, as listed in Table 2-9. 
• 
Wave data for' all the ,known tests in water 
with HoE. charges larger than 10 lb, as list- 
ed and referenced in Table 2-100 
Surface. waves 
produced 
by 
an 
underwater 
explosion are generated by pulsations of the 
cavity fonned at the surface when the oscillating 
, :,..;~:::~ '=::.:!:: ~~~ough the surface (vents). 


• 
At a short range, within a few cavity 


ra~the fLrst wave crest is so high and steep 
that it is unstable', and it spills forward 
turbu~ 


lently and dissipates much of its energy. How- 
ever, as the wave train expands and attenuates, 
the waves become shallow and smooth, their 
number increases, their energy remains nearly 
constant, and the highest wave appears succes- 


si~ely later:in the train. A wave gage record for 
such an expanded wave train is shown in Figure 
2-106. The height H ofthe highest wave (trough· 
to-crest) de,creases nearly inversely with range R ; 
but the length, speed, and period all remain 


nearly constant for. this peak wave, which is the 
one temporarily located at the flIst anti-node of 
the wave envelope as the waves move back in the 
envelope. 


_ 
An 
approximate 
mathematical theory 


h~n developed for shallow waves propagat- 
ing from a shallow disturbance on the surface of 
water with constant depth. The results of this 
theory compare well with measured wave trains 
when a motionless "pseudo-<:avity'; is selected to 
approximate the wave generator or initial condi- 
tions. The "pseudo-cavity" is a shallow motion· 
less depression of the water surface near the 
explosion. Its shape is derived empirically to cor- 
respond with the waves measured far from an 
explosion, but it does not necessarily correspond 
accurately with the actual shape of the water 
~e at any time. 
. 
, 


• 
All of the known measurements of peak 
wave height over deep water for charges greater 
than 10 Ib are plotted in Figure 2-107.' The verti- 
cal bars indicate the variation of HR for repeat- 
ed shots and/or various gage ranges Ro There is 
,little discernible effect of depth of burst except 
near the "upper critical depth," where the H.E. 
charges 
were ·only 
partially 
submerged 
and 
where the peak wave height nearly doubled. Al- 
though current studies indicate that the upper 
critical depth effect may exist for nuclear as well 
as for H.E. charges, there is no fmn evidence to 
support the inference as yet. 
II A plot comparable with Figure 2-107 of 
the data for H.E. charges smaller than 10 lb 
exhibits a pronounced upward hump in wave 
·height versus depth of burst at~ lit: 3~{ (and 
possibly also at db 
Ill:: 8k1ii). This increase in 
wave generation efficiency at the "lower critical 
depth" occurs when the expanding underwater 
cavity ruptures the layer of water above it (or 
vents) just before it becomes fully expanded. 
For large submerged explosions, the expanded 


, bubble becomes tangent to the surface when the 


~ 


Table 2-9. 
_ 
Measured Wave Data from Nuclear Tests II 


NUMBER 
WATER 
DEPTH 
PEAK WAVE- 
CAVITY 
YIELD 
DEPTH 
- OF 
HR/2 x 104.JW 
SERIES 
AT CHARGE BURST HEIGHT x RANGE RADIUS dw /l00(W)1/4 
,SHOT, 
W (kt) 
" HR/2,000 (ft)t 
dw (ft) 
~ (ft) 
Rc (ft) 
min - max 


'7 JJP 


100(W)1/4 
miD- max 


! )-:lL1) 
I • 


140 
140 
• • - 


\..' 
HARDTACK 
UMBRELLA 


2 .i 
3,006 
500 - - • - 
HIRDTACK 
WAHOO· 


3 
23.5 
180 
90 
36 - 60 
1,000:1: 
0.82 
0.74 - 1.24 
CROSSROADS 
220 
Shallow 
BAKER 


4 
32 
15,000 
2,000 
118 
490t 
63 
2.08 
WlGWAM 
238 
Very Deep 


5 .i 


120 
0 
• • .- 
REDWlNG 
FLATHEAD 


6 • 


140 
0 


.- 
• • - 


,. REDWlNG 
DAKOTA 


7 • 


230 
0 
• ..- 
REDWING 
NAVAHO 


8 
.7,000 
145 
7 
. 232 
- 312 
I,OOOt 
0.16 
0.16 - 0.28 
CASTLE 
916 
Shallow 
UNION 


9 
13,500 
220 
7 
.426' . - 438 
I,OOOt 
0.20 
0.37 - 0.38 
CASTLE 
1,080 
Shallow 
YANKEE 


.-H is twice tbe measured beight of the' peak crest except for shots 8 and 9 where H is the measured height 
of the first crest from the foUowing trough. H is corrected for uniform water depth = db by Green's law. 
tValue deduced· from measured surface wave train (Kaplan and Wallace, see bibliography) is a lower limit 
considerably sm&uer than actual values, which are unknown. 
tMeasured values of the column radius (Young, DASA 1240-1(9», see bibliosraphy. 


Table 2-10. '_ Measured Wave Data 'from Large H.E'. Tests II, 


min - max 
W (kt) 


YJELD 
WATER 
DEPTIl OF PEAK WAVE CAVIn' 
HR;2 x 1O"V'W 
t- 
-+-D_EPTH_'-+--:-B_URS_T_+-_HE_IG_HT--:-4 
RA-...:.D_IU 
--lS d../l00(W)1I4 I- 
~ 
d"" (ft) 
~ (ft) 
HR/2 (ft)2 
.Re (ft) 
SHOT 


125 Ib 


385 Ib 


WES Test Series at Lake Ouachita (from Pinkston and from Whalin) 


0.625 x 10-4 
100 
0 - ,25 
80 - 
110 
22.5' 


1.97' ~ 10-4 
100' 
0 - 36 
147 - 
225 
17.0 


.90 - 1.24 


1.25 - '1.91 


WESTeSt Series at Mono Lake (1965 - 1966) 


(prellininary data by private cOmmunication with R. Whalin and M. Pinkston and from Walter) 


I 
1965 - 6 
1 0.46 
x 10-2 
130 
0 
1650 - 1900 
68 
5.2 
2.4 
- 2.8 
(9,200 Ib) 


-2 
.67 
2100,., 3900 
69 
3.1 
- 5.7 


- 4 
1.04 
2100 - 2850 
71 
3.1 
- 4.2 


• 


- 3 
1.40 
1750- 3000 
74 
2.6 
- 4.4 
() 
-5,7,8,9 
9 - 51 
61 - 80 


j 
1966 - 3 
=-130 
5.2 
1150 - 1250 
1.7 
- 1.8 


I 
- 2 
! 


23 
850 - 1200 
1.2 
- 1.8 


n 
- 1 
60 
1300 - 1400 
1.9 
- 2.1 


ij 
HYDRA Ita (Van Dorn) W = 10" Ib HBX 


H 


- 7 : 0.723 x 10-2 
300 
140 
1650 - 1970 
10.3 
1.94 - 2.31 


Ii 


(14,450 Ib) 


!t 
- 8 . 
80 
,975 - 1340 
U5 - U5 


~ 


-9 i 
35.5 
1260 - 1590 
1.48 - 1.87 


- 10 
7.1 
1610 
1.88 


t' 
-11 
15.9 
1510 - 1840 
1.77 - 2.16 


- 12 
159 
1490 
, 
, 1.75 
- 
- 13 
17.7 
1450 
1.71 


WATER LEVEL 
L 


Figure 2- 106.II Wave Gage Record for an Explosion- 
Generated Wave TrainII 
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• 
burst 
depth 
db 
equals 240W' /4. The lower 
critical depth is somewhat smaller at about db = 
170W' /4: However, the effect on wave genera- 
tion is much smaller for the larger charges in 
Figure 2-107, and the height of the peak wave 
generated by large explosions in "deep" water 
can be approximated simply by 


HR/W~i4 
= 16 ± 4, 


where H is the amplitude of the maximum of 
the first envelope trough-to-crest in feet,· R is the 
radial distance from the explosion in feet, and 
WH E is the charge weight in equivalent pounds 
of TNT. An equivalent approximation is 


-HR = 40,500 (W)0.54 
± 25 percent, 


where W is the yield in kilotons and the other 
symbols are as defined above. 
Similarly, the 


. length and period of the peak wave can be ap- 
proximated by 


L = 15.5 
(WHE)O,288 ft, 


T = 1.74 (WHE ) 0.144 
sec. 


. These three relationships were derived empirical- 
ly from 
test data for 
125- and 385-lb H.E. 
charges and the two nuclear explosions shown in 
Figure 
2-107. For a pseud<Xavity ofa given 
shape, 
which 
has 
a 
potential 
(gravitational) 
energy comprising a fixed fraction'of the explo- 
sive yield, all the dimensions of the cavity and 
the wave train are proportional to W:/~ or Wi /4 
(Froude scaling), and the exponents in the three 
equations given above'are 1/2, 1/4; and 1/8, re- 
spectively. Since the empirical exponent in the 
flIStequation was found to· be 0.62 for small 
H.E. charges, it Seems reasonable to suppose that 
the "theoreticaJ" value of 1/2 may be the best 


value to use for extrapolation to large yields (the 
high kiloton and megaton range). With this ex· 
ponent, the data in Figure 2-107 have been re- 
plotted as a function of reduced water depth 
(rather than charge depth) for use in Probiem 
2-35 (Figure 2-128). 


_ 
The horizontal line' in Figure 2-1 28 can 
be used to predict the peak wave height in deep 
water. However, it can be derived by the follow- 
ing simple procedure to illustrate qualitatively 
the 
physical 
process of wave 
generation 
in 
"deep" water to define "deep." The approxima· 
tion is made that an explosion in deep water at 
any 
depth of bUr.'t above the lower critical 
depth generates a parabolic (nearly hemispheri- 
cal) cavity upon venting, and that the potential 
energy of this cavity Ec\$equal to 5 percent of 
the explosive yield Eo' That is: 


Ec = ; Pw gR: = 0.05Eo 
= 0.05 x 2.86 x 1012 W ft-lb. 


_ 
Thi~ gives the following cavity radius, 
which is taken equal to the cavity depth and the 
radius of the plume (or column) of water pro- 
jected upward: 


R c = 256 (W)I/4 
ft. 


and if the density of seawater and the hydrody- 
namic yield of the explosion are taken as 


P9l g = 65 ~ 
and e 
= 1,025 cal. 
ft 3 " 
0 
gm 


This cavity radius is not much larger than the 
bubble radius and depth of burst corresponding 
to-the 10wercritica1 depth. The restriction on 


1 
! 


• 
water dl'pthdw in the equation for Rc is taken 
as the minimum water depth that can be con- 
sidered "deep" i.e., the depth for which the 
assumed cavity shape touches bottom as it vents. 


_ 
The dissipative effects close to the cavity 
are approximated by the assumption that the 
leading peak wave crest spills forward and breaks 
turbulently at R =2Rc to the limiting height of 
a wave that can propagate in a stable manner 
over deep water (the MicheU limit, or Hmall = 
0.14L). FinlWy it is assumed that the amplitude 
'of the wave envelope H is attenuated inversely 
with range R as it continues to expand without 
dissipation. 
This gives the 
"prf'dicted" !,eak 


wave height, when the water depth at the charge 
location is "deep" (d", ';;> 256WI / 4 ) as 


over shallow water (Table 2-9) and which will be 
discussed iri succeeding paragraphs. 


R = 100 (W)1/4 if d' 
<; R 
e w e 


Once again, the dissipative effects Close' to the 
cavity are approximated by the assumption that 
the leading peak wave breaks down to the limit- 
ing height of a stable progressive wave, but, for 
shallow water, the approximation is that Hm all 
= (3/4) d", = Miche limit at R = 2R c ' If this peak 
wave height is then attenuated inversely with 
range R, the "predicted" peak wave height when 
the water depth at the charge location is "shal· 
low," (dw <;; 100WI /4 ), becomes 


HR = 49,000 ~ 
HR = ISO d 
Wl/4 


'" 
' 


o 


which gives the line with the 4S degree slope in 
Figure 2-128. This "prediction" can be used for 
any depth of burst (0 <; db <; d", las shown by 
measured data for H.E. charges in shallow water. 
As shown in Figure 2-128, there is much scatter 
in the nuclear wave data in shallow water and 
some of the measured peak wave heights are 
more than twice as large as "predicted" by the 
foregoing equation, which should be regarded 
, only as a rough rule-of-thumb for correlating the 
measured data in Figure 2-128. 
' 
__ The only actual measured data for 
Hie radii of the water columns (or cavities) for 
nuclear ,?ursts near the surface of water are for 
shots BAKER and UMBRELLA (Table 2-9). 
These values of Rc and those measured for sub- 
merged H.E. explosives are roughly equal to the 
value given by the equation for maximum bub- 
ble radius in para 
i 
h 2-67, which gives a burst 
, estimate for R 


which corresponds to the horizontal,line in Fig· 
ure 2-128. This "predicted" value 
~n be used 
for depths of burst ranging from zero to at least 
five 
times the 
lower critical depth (at d 
== 


170W1 / 4 ) judging by the comparison with WIG- 
WAM data, ~ven though the cavity radius becomes 
too large for 
depth~ of burst much below the 
lower critical depth. 
, 
11 ~e ~e ~ Figure 2-128 can be e~tend­ 


"'.......u 
~iGllUW 
water by a comparable sunple 
procedure. ,In shallow water the cavity both 
touches bottom and vents to the atmosphere 
early during its expansion. Hence, much less 
work is done on the water by the expanding 
gaseous 
bubble, 
and 
much 
less 
energy 
is 
propagated as water surface waves over shallow 
water. Experimental data indicate that explo- 
Sions in shallow water generate a dry cavity sur- 
rounded by a nearly vertical wall of water which 
extends upward as a thin lip to form 
the 
"column.'" A lower limit 
for the maximum 
radius of this cavity (or column) is given by the 
{oUowing relationship, which was derived from 
the t.nergycalculated from measured wave trains 


~~o 
2-'228 


, 
iJ 


, 
f1 


l, ~. 
• 


n so , 
c argesdetonated at the 
s 
ace generate crater radii approximately twice 
as large as nuclear explosions of the same yield, 
so the nuclear-TNT equivalence is only about 10 
percent. The column radii observed for TNT det- 
onations near the surface of the water are rough- 
ly equal to 470W1 / 3 , ,which is the bubble radius 
given by the equation in paragraph 2-67 for a 
surface burst. This occurs even when the bubble 
radii exceed the water depth by a factor of 5, 
Therefore, it might be expected that a nuclear 
burst at the surface of water would produce a 
column radius about half as large as the H.E. 
value, or 240W1 / 3 . This latter value is much 
larger than the lower limit estimated above (R c 
=IOOwl /4). This lower limit is also exceeded by 
the crater radii measured in the bottom after 
explosions in 
shallow water, 
and 
the 
water 
column radii probably exceed the crater radii. 
Although there are no substantiating data, nu- 
clear explosions at the surface of water of any 
depth probably generate water columns with 
radii between lOOW l / 4 and 240W l / 3 feet. 


2-78 
Refraction and Shoaling II 


• 
Explosion-generated waves, in common 


wilr'"wind 
waves 
and 
tsunami (earthquake- 
generated) waves, are affected by changes in the 
depth over which they propagate. As the waves 
move into water shallower than about one-third 
their length 
(or shoal), their period remains 
fixed, but both their speed and their length de- 
crease 
(between 
successive crests) and 
their 


height first decreases about 10 percent and then 
begins to increase'. Because of the change in 
speed, a wave bends or refracts if it moves into 
shallow water at an angle to the bottom con- 
tours. Refraction can increase or decrease the 
localized wave height, depending on the hydro- 
graphy over which the wave moves. No general ' 


methods for determining refraction effects are 
given here; however, such techniques for calcu- 
lating these effects are available ("Oceanograph- 
ica1Engineering," 
R. 
L. Wiegel, see bibliog- 
ra. 
' 
The increase in height of a wave relative 
to Its length (steepening) continues until the 
wave becomes unstable and breaks, unless the 
bottom slope is so shallow that bottom friction 
dissipates the wave before it breaks. In shallow 
water without dissipation, steepening increases 
the wave height as the inverse fourth root of the 
water depth. Many types of waves become un- 
stable and break (or spill) when their height 
exceeds either 75 percent of the water df'pth 
(the Miche limit) or 14 percent of the wave 
length over deep water (the Michell limit) as de- 
scribed by Weigel (see bibliography),' 
'I UNDERWATER CRATERING II 


_ 
If an explosion occurs in or even close to 
a layer of water overlaying bottom material, a 
significant crater forms in the bottom material 
whenever the gaseous bubble or water cavity 
formed by the explosion contacts the bottom. 
An underwater crater is similar to. an onshore 
crater forined by an explosion near the ground 
since both are characterized by a dish-shaped 
depression, wider than it is deep, and surround- 
ed by a lip raised above the undisturbed surface. 
In the case of most underwater craters, however, 
the observed ratio of crater radius to depth is 
larger and the lip height is smaller than for cra- 
ters from comparable bursts in similar bottom 
materials onshore. These differences are caused 
by water displaced by the explosion washing 
back over the crater. This flow increases the cra- 
ter radius by as much as 10 percent and de- 
creases the crater depth by as much as 30 per- 
cent. An exception to this general rule occurs 
when the water layer is so shallow that the lip 
formed by the initial crateringextends above the 
surface of the water. Such craters (termed "un- 


.n'"1 
Z.aL 


I 


• 


II 
. 


washed craters") approach onshore craters in 
appearance, with higher lips and larger depth-to- 
radius ratios than washed craters. This section 
describes only the crater characteristics directly 
related to the explosion. The effects of tidal cur- 
rents and the collapse of unwashed crater lips by 
hydrostatic forces are not considered. 
"The characteristics of underwater craters 


=;:~ ;!'e~icted at least qualitatively from the . 
depth of burst, water depth, bottom composi- 
tion, and weapon yield. The use of empirical 
curves reIRting the depth and radius of the crater 
for surface and bottom bursts over a clayey-sand 
bottom is described in Problem 2-36. Predictions 
for other bOttom materials can be made by ap- 


plying tabulated factors to the values predieted 
for clay-ey sand, as described in Problem 2-36. It 


. is evident that larger craters are produced by 
moving the charge toward the bottom; increas- 
ing the water depth above a bottom burst (up to 
a point); or decreasing the water depth beneath 
a surface burst. 
_ 
A very large yield is required.to form an 


un:ed crater in water having a depth that is 
common to most harbors. For example, with 30 
£t of water over clayey-sand, 4 Mt would be re- 
quired on the bottom to produce a dry crater 
lip. It is quite unlikely that materials. that tend 
to. flow, such as oceanic ooze or fme-grained silt, 
could support the hydrostlltic pressure exerted 
on a dry crater lip. 


"f, 


Problem 2-28. 
Calculations of Underwater 


Explosion Phenomena 
II Calculations .can be performed to deter- 
mine 
bubble characteristics described in para- 
graph 2-67. The characteristics include the maxi- 
mum radii, bubble periods, and migration be- 
tween 
maxima for. events to the end of the 
second bubble period. The calculation methods 
do not apply to cases where the bubble touches 
the bottom, and the methods contain uncertain- 
ties not yet fully 
resolved when the bubble 
center is very close to or above the surface. Fur- 
thermore. predictions that are made for events 
increasingly 
removed 
in 
time 
from 
the first 
maximum 
have 
increasing 
uncertainty. 
Two 
methods for calculating these bubble character- 
isLcs are demonstrated in the example. The first 
method uses the equations given in paragraph 
. 2-67 and Figures 2-108 through 2-110. The sec- 
ond method uses equations introduced in the 
example and Figure 2-11 I. The first method is 
graphical and is less accurate than the second, 
bu t it provides a quick estimate. 
II Example 
- 
. 
' 
Given: 32 kt= at a depth of 2,000 ft in 
very deep water. 
Find: 
a. 
The initial bubble period and radius, 
b. 
The 
depth 
of the second maximum 
radius, 
c. 
The second bubble period, 
d. 
The radius at. second maximum, 
e. 
The depth of third maximum· radius, 
f. 
The third bubble period and radius at 
third maximum. 
Solution (Method-}): 
a. 
Calculate the initial radius A I and the 
period TI , by using equations given in paragraph 
2-67: . 


wl/3 
Al :: J'-- 
ZI/3 ' 


Al 
:: (1,500) ( 
(32)1/3 
) 
:: 376 ft, 
(2,000 + 33 )1/3 


Wl/3 
TI 
:: K'-- 
Z 516 ' 


TI 


:: (SIS) ((32)./3 
) 
:: 2.86 sec. 


. 
(2,033)516 


b. 
Find the reduced hydrostatic head: 


ZI 
. 2,033 
- 
::--= 5.41. 
Al 
376 


. From Figure 2-108, for n :: I, 


.:121 
-- = 0.28, 
ZI 


.:121 :: (0.28)(2,033) = 569 ft 
(migration during fIrst period), 


d 2 = 2,000 - 569 :: 1,431 ft.· 


., Note that rellllts are ellSitiwi to IBPh readinls. 


. i 


c. 
From Figure 2-109, for n =2, 


T2 = 0.90 T1• 


T2 = '(0.90)(2.86) = 2.57 sec. 


d. 
From Figure 2-110, for n = 2, 


A 2 = O.i17 Al = (0.77)(376) = 290 ft. 


, 
e. 
From Figure 2-108, using Zl!A 1 =5.41· 
and n =2, 
. 
a. 
By the same procedures used in Method 


~2 
-- = 0.23, 
ZI 


I, 


Ai = .376 ft, 


liZ2 = 0.23 ZI = (0.23)(2,033) = 468 ft, 
'T1 - 
2.86 sec. 


b. 
Using the first reduced equation given 
above, 


d 3 = 1.4,31 - 468 = 963 f1. 


f. 
From Figure 2-109, for n = 3, 


~1 
= 0.277 ZI' 


~1 
= (0.277)(2,033) = 564, 


( ~l) = 
. 
(~)3f2 
ZI' 
3.5 
2,033 
T3 := (0.66)(2.86) = 1.89 sec. 


From Figure 2-110, for n =3, 
I 


A 3 = 0.47 Al = (0.47)(376) = 177 ft. 


I 
d 2 
=2,000 - 564 = 1,436 ft. 


Solution (Method 2): This method uses 
Figure 2-111 to obtain values for bubble energy 
ratio, 'a/'a.il' and the following reduced equa- 
tiou::.; 
, 


c... The second bubble period T2 can be 
obtained from the third reduced equation given 
above, and T) from (a) after fmding the energy 


tu2 = 0.303 Z2 = (0.303)(1,469) = 445 ft, 


II 
. ratio between the flfst and second bubbles. Us- 
ing the reduced migration tul/Z I = 0.277 to 
en ter Figure 2-1 I I, the energy. ratio is , 2!'I = 


, 0.324. Thus, 


_( T2 ) = (2,033~/6 (0.324)1/3, 
T1 
1,469) 


(~2) = 


(~')= 


3 5 (~)3/2 (0324)112 ( 2,033)2 
. 
2,033 
. 
1,469 
' 


0.303, 


T2 = (090)(2.86) = 2.57 sec. 


d. 
The radius of the second maximum is 
obtained 
from 
the 
second· reduced equation 
given above. 


(~).(~t' (0.324)1/3, 
Al 
1,469 


(~:) 0:765, 


A2 = 0.765 Al 
= (0.765)(376) = 288 ft. 


e. 
To 
calculate the depth of the third 
maximum, it is flfSt necessary to find the reduc- 
ed migration tu2/Z 2 . From the flfSt reduced 
equation, 


d 3 = 1,436 - 445 = 991 ft. . 


f. 
From c., '2/'1 = 0.324. Using the re- 
duced migration tu2/Z 2 = 0.303 to enter Fig- 
ure 2-111, the energy ratio '3/'2 =0.158. 
. 


( :: ) = (0.324)(0.158) = 0.05 f2. 


From the third reduced equation 


(_ 
T3 ) 
__ (2.033)516 
- 
(0.0512)1/3 = 0.658, 
T1 
1,024 


T3 = (0.658)(2.86) = 1.88 sec. 


c. The second 
T2 =2.57 sec 
2.57 sec 
bubble period 


d. The radius at 
A 2 = 29U ft 
288 ft 
second maximum 


e. The depth of 
d 3 = 963 ft, 
991 ft 
third maximum 
radius 


II 
From the second reduced equation, 


, 
I 


( A 3') = (2,033)1/3 
(0.0512)1/3 = 0.467 
AI, 
i 1,024 
. 


A 3 =(0.467)(376) = 176 f1. 


'f. The third 
bubble period 
and radius at 
third maximum 


T3 =1.89 sec 
1.88 sec 
A 3 = 177ft 
176ft 


.Merhnd J 
Method 2 
Al = 376 ft 
376 ft 


.T 1 =.2.86 sec 
2.86 sec 


Answer: 
a. The initial 
bubble, period 
and radius 


b. The depth of 
d 2 = 1,431 ft 
1,436 ft 
the second 
maximum radius· 


_ 
Reliilbiliry: 
The 
uncertainty 
in 
the 


va~for maximum radii calculated by the 
methods described above are about +20 percent 
and -35 percent, based on very limited informa- 
tion. 
II Related Materiill: See paragraph 2-67. 
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Figure 2-108. II 
Reduced Nuclear Bubble Migration 
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Successive Nuclear Bubble Period Ratios • 
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Figure 2-110. 
• 
Successive Nuclear Maximum BUb~le Rad,ius Ratios II 


~) 


r---;...-.--..,...------r-----~r------ ~ 


It') 
0 


('oj 


II 
M 
c: 
II 
- 


c: 


ci 


III0 


.~ 
ci 
C 


.<11 


~ 
N 
a:: 
> 
~ 


Z· 
II> 
C 


0 
w 


i= 
II> 
M 
« 
:c 
0 
a:: 
J:J 
!:? 
:l 
lI:I 
:e 
.. 
<II 


0 
.11> 
.r 
o 
, 
w 
~ 
. 
) 
U 
:::l 
:> 
z 
'"" 
0 II 


('oj 
w 
d 
a:: 


1 
..:. 


I 
...... 
~ 


II>.. 
5. 
0 
ir· 


L-__"'--__.1-__.1-__.L-__...L-__:3:::...._.J..__.J 0 


~ 
.~ 
~ 
c; 
0 


Problem 2-29. 
Calculation of Shock Wave Parameters 
for Free Water (Deep Bursts) 
II Figures 2-112 through 2-116 sh'ow the 
parameters for the primary shock wave from an 
underwater explosion in free. isovelocity water 
(no reflections, no refractions). Figure 2·112 
shows the peak pressure as a function of slant 
range for various yields, where the peak pressure 
is lower than 3.000 psi, and as a function of 
reduced slant range (R/W I/3 ) for peak pressures 
above 3.000 psi. Figure 2-113 shows the shock 
wave impulse as a function of sIan t range for 
various yields. while Figure 
2-114 shows the 
shock wave energy flux as a function of slant 
range for the same yields. These latter two fam- 
ilies of curves show values of impulse and energy 
!lux integrated to a time of 6.78. where 8 is the 
time constant for the shock wave. which is de- 
fined as the time for the shock wave pressure to. 
fall 
to approximately 
37 percent of its peak 
value. 
_ 
The time constant 8 is shown in Figure 
2- ~as a function of slant range for the same 
family 
of yields shown in Figures 2-113 and 


2-114. 
Figure 
2-1 16 
shows 
a dimensionless 
pressure-time curve. Using values of 8 from Fig- 
ure 2-115 and values of Pm from Figure 2-112, 
Figure 
2-1 16 may 
be used 
to construct the 
pressure-time shape of the shock wave in free 
water at various ranges from underwater explo- 
siWf various yields. 
. 


• 
Scaling. 
For. yields, other 
than trose 


shown in Figure 2-1 12 through 2-1 15. use linear 
in"ero:ation between appropriate curves. 
Example _ 


Il'en: A 50 ~urst at a depth of 1,000 ft 
in deep water. 
Find, 
The peak shock pressure Pm' the 
shock wave impulse I, energy flux E. and time 
constant 8 at a 1,000 ft depth. 4,000 yd from 


the burst. 
Solution.: The desired shock wave param- 
eters can be read directly from Figures 2-1 1:2 
through 2-1 15. 
Answer: 


Pm 
:::: 470 psi, 


1 
=:::: 28 psi-sec. 


E 
=:::: 
103 in.-lb/in.2 , 


8 
=:::: 50 msec'. 
III Reliability' 
Values 
shown 
in 
Figures 


2" I 12 through 2-1 15 are estimated to be accu- 
rate within ±10 percent for shock pressure. and 
±20 percent for the other parameters for yields 
between I kt and I Mt. For yields below I kt or 
above I Mt, the values are more uncertain. The 
assumption 
of isovelocity, water becomes in- 
creasingly uncertain for ranges beyond a few 
k-,'IOards in the ocean. 
Related Material: 
. 


ee paragraphs 2-68 through 2-72. 
Figures 2-112 through 2·1 15 are based on 
the following equations: 


. 
(W1/3)0.91 
I (to 6.78) =1.176 x 104 wl/ 3 R 
psi-sec! 


E(to 6.78) =3.976 x 109 WI/3(WI/3)2.04 in.-Ib 


. 
. 
R 
in.2 


2-"283. 
~/ 


The values of these 
parameters are re- 
lated to those, for TNT (Cole, 1948) by an ap- 


II 


( 


113)-0.22 
e=.2.2"/4 wI! 3: 
k~ 


2-~ 
~z. 


msec.' 


proximate equivalence factor. Snay and Butler 
(1957) (see bibliography) report that an under- 
water nuclear burst with a radiochemical yield 
of 1 kt gives approximately the same shock wave 
parameters at the same distance as a TNT charge 
weighing 0.67 kiloton (where 1 ton =2,000 Ibs). 
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Figure 2-112... Shock Wave Peak Pressure in Free Isovelocily Water. 
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Figure 2-113. II 
Shock Wave Impulse vs Slant Range in 
Free lsovelocity Water • 
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_ 
Shock Wave Energy Flux 
V$ Slant Range 
,in Free Isovelocity Water 
• 


SLANT RANGE (kml 


~ 
i 


i 
I 


.J 


I 


10 


SLANT RANGE 
hhousands of Vds) 


£'~. 
~ ... 


1 
,0 1 
SLANT ,RANGE (km) 
10 
100 


100 
10 


SLANT RANGE Cthou.ndI of ydsl 


2-115. _ 
Shock Wave Time Constant in Free lsovelocity 
Water vs S1lnt Ringe • 
Figure 
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. Figure 2-116. _ 
Pressure-Time fo~ Free lsovelocity Water 


. Shock Waves .fromNuclear Expl.osion 
_ 


Problem 2-30. 
I 
Calculation of Shock Wave Parameters for Reduced Water 


. Depths, ~wlW1/3 <; 2,000 ft/kt1/3 


• 
Figu'res 
2~ 117 through 2-120 show the 
influence 
o~ surface and bottom reflection on 
the primary shock wave. The curves are drawn 
through values of peak pressure and impulse at 
mid-depth 
both 
for explosions at mid-depth 
(Figures 2-117. and 2-118) and at the bottom 
(Figures 2-1:19 and 2.'120). The solid curves indi- 
cate the approximate coverage of the data. As 
the water depth approaches 7.ero the dashed 
curves 
approach 
the 
values 
for 
air blast. at 
ground level. Transition to .the free-water shock 
__.parameters 
(Figures 
2-112 
through 
2-116) 
occurs 
as 
the 
water depth increases. In the 
anomalous region near the water surface, the re- 
flected tension waves lead to a sharp reduction 
Ji! 


f 
re~sure and impulse (Figure 2-94). 
The curves in Figures 2-117 through 
2- 
ere derived mainly from H.E. data; how- 
ever, they were shifted to agree with nuclear 
data for th'e limit of free isovelocity water. In 
most· instances, the shift (or conversion factor) 
was in agreement with the factors determined 
separately for converting from H.E. to nuclear 
rl;\t::t in deeo water. Most of the data were ob- 
tained. for the mid-depth location of both the 
charge and' the pressure gauges. All of the curves 
pertain to either a coral sand or an ordinary san,d 
bottom material. 


shave 
een made over concrete, sand, and clayey-silt 
bottoms to assess the effects of shock reflection 
I 
fr0m variou!' bottom materials. Measurements 
were made within 48 ft from 1Golb TNT charges 


~8 
2-24Q. 


in 20 ft of water. Both the charge and gauges 
were placed deeper than mid-depth to minimize 


• 


effects of surface reflection. 
Briefly, the. following effects were ab- 


. 
e 
or the shock strength, impulse, and ener- 
gy flux as compared with values at the same 
range in free water. Near the bottom and away 
from the charge, the pressure fronts of the inci- 
dent and reflected shocks were nearly equal. 
Near the charge and away from the bottom, the 


renect~d shock was about one-fifth as strong as 
the incident shock. The impulse of the shock 
wave increased as much as 100'percent near the 
charge over concrete, but there was little change 
at the greatest range from the charge. The im- 
pulse increased near the charge .but decreased 
further away from it overclayey-silt and sand. 
The increase was as much as 40 percent over silt 
and 20 percent over sand. The decrease was as 
much as 50 percent over silt and 30 percent over 
sand. The observed changes in energy flux were 
roughly comparable with those for impulse, al- 
though the energy decreased as much as 30 per- 


~om the charge over concrete. 
~ The shock wave duration T in milli- 
seconds can be approximated with 


KJ 
T=- 
Pm 


K l'l:: 1.8 + 0.0021 
(dVo' /Wl(3) 


for 2,000 > (dw /WI'3) > 100 


K l'l:: 3.0 - 0.01 (dw/Wl/3) 


Use air blast loading 
a 


Example II 
, 
Given: - A 
1- Mt explosion at a deptt 
2,500 ft in water having a depth of 5,000 f1. 
Find: The peak pressure, the impulse. 
shock wave duration,and the energy flux 
depth of 2,500 ft and a horizontal rangf 
20,000 ft. 
Solution: The reduced water depth is 


I = shock wave impulse (psi-msec) 


where 


P' =peak shock pressure (psi). 
m 


dw = water depth (ft) 


The shock wave energy flux for an 
explosIOn at a reduced depth of approximately 
250 ft(kn l ,' 3. is related to the shock wave im- 
pulse I, and peak pressure by 


dw 
5,000 
- 
= --.,;~:..::-- = 500 ft- 
WI/3 
(] ,000)1/3 


The reduce.d burst depth is 


\ 


and at reduced depths equal to or greater than 
:2.000 ft(kt)113 by 


2,500 
(1,000)1/3 


The reduced range is 


= 250 ft. 


R 
20.000 


W1/3 = (1,000)113 = 1.000 f1. 


For intermediate depths, a linear interpolation 
between the two values can be adopted: 


_ 
The following approximations .can 


~stimate the shock parameters at van· 


ous gauge depths when -the charge is at mid- 
. depth. when 250 < R/W113 <.2,500 and when 
70 <d/W' /3 < 200. 
w 


Since the explosion is at mid-depth,' Fig 
2-1.17 and 2-1 I 8 apply. From Figure 2-117, 
peak pressure is 


Pm = 700 psi. 


From Figure 2-1 18, the reduced positive imp 
is 


'Since 2,000 > (dw /Wl/3)> lOa, the shock ~ 
duration equations give 
. 


K = 1.8 + 0.0021 (dw /WI/3) 


Gauge Depth 
'Fraction of Value 
Water Depth 
at Mid-Depth 


f9 


for peak pressure 


3/4 to .1 
0.8 for impulse 


0.7 for energy impulse 
fa 


for peak pressure 


1/4 
0.7 for impulse 


0.6 for energy nux 


I 


W1/3 = 9,000 psi-msec/(kO'/3. 


~ 


K = '1.8 + (0.0021 )(500) = 2.85 


Answer: The peak pressure is 


i 


The shock wave energy flux is 


I = (9,OOO)(W I/3 ) = 90,000 psi-msec 


Pm = 700 psi. 


The positiveimpLilse is 


The shock wave duration is 


T=KI 
= (2.85)(9 x 104)= 


Pm 
700 
366 msec. 


£250 = (1.2 
X 1()"'4 )(700)(9 x )Q4) 


= 7.6 X )03 in. Ib/in.2 • 


...... Reliability: 
Peak 
pressures 
obtained 


~igUres 2-) 17 and 2-) 19 are estimated to 
be reliable within about ±30 percent. Impulses 
obtained from Figures 2-) 18 and 2-) 20,and en- 
ergy fluxes are estimated to be reliable to within 


. a~. 
ut~±50 percent. 
' 
~ Related Material: 
See paragraphs 2-68 
through 2-72. 
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Reduced Positive Impulse of Shock Wave vs Reduced Water 
Depth. 
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. Problem 2-31. 
Calculation of· Peak Overpressure Along the Water Surface 
as a Function of Underwater Burst Depth 


_ 
Fi2ure 2-121 shows the relations to ob- 


tai~ max·lmum air blast overpressure from a 
scaled underwater explosion as a function of 
depth of burst and horizontal distance. The air 
blast is for an altitude just above the surface of 
the water. To detennine the air blast at other 
altitudes. see the references listed in paragraph' 
"-7'"' 
. 
. 
___ 
E.~ample _ 
~. A 10 kt explo~ion at a depth of 500 
ft in deep water. 
Find: 
a. 
The horizontal distance at which a 0.5 
psi peak air blast overpressure would be pre- 
dicted. 
b. 
The upper and lower limits of the range 
for item a. that would cover the estimated accu- 
racy. 
Solution' The reduced depth of burst is 


500 
= 
= '232 ft/(kt)I/J. 


(10)1/3 
/ 


From Figure 2-121, the intersection of the 0.5 
psi curve with a reduced depth of 232 ft/(kt)1/3 
occurs at a reduced horizontal distance of 0.7 
kft /(k t)113 . 


Answer:' 
a. 
The 
predicted 
horizontal distance at 
which 0.5 psi will occur is 


. d = 0.7 
Wl/3 
= (0.9)(10)1/3 


d = 1.50 kilofeet = 1.500 feet. 


b. 
The estimated accuracy in predicting 
the distance for a given overpressure is ±30 per- 
cent (see "Reliability" below). Thus, the limits 
are 


. d = 1,500 x 1.3 = 1.950 ft (upper limit). 


d = 1,500 x 0.7 = 1,050 ft (lower limit). 


_ 
Reliability: The accuracy of the relation 
in Figure 2-121 has not been tested rigorously. 
The relations provide an approximation based 
on high-explosive underwater tests and canfinn- 
ed by a few measurements made at nuclear 
underwater tests. However, it is felt that ±30· 
,percent in predicting the distances for maximum 
overpressure is a realistic accuracy estimate. 
• 
Related Material: See paragraph 2-72.. 
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Problem 2·32. 
Calculation of the Extent, Initial Vertical Velocity, 
and Approximate Height of the Spray Dome 


_The maximum extent of the spray 
~e estimated ",:ith the followiryg rela- 
tions. 
which are valid 
for reduced depths of 
burst as large as: 


lated from: 
Initial spray velocities can' be calcu- 


Pm 
Yo = (144 cos' 0) -U 


o 
= 85.3 - 0.021 '(db /W I/3 ) 


ma~ 


where (as shown in the following sketch) 


where at a point on the surface. 


Pm = peak shock pressure (psi). 


U = shock front propagation velocity. (ft/sec). 


Yo 
= initial spray velocity (ft/sec.) 


hi = the weapon yield (kt). 


r,d 
= ,the horizontal extent of the spray dome_ Values of Pm are given iri. Figure 2-112, 
(fll. 
Values of U are given in Figure 2-112 as a func- 
tion of Pm' 
In practice. the presence of the 
o= the spray dome 'angle. the angle between 
bottom does not appear to have an influence 
a vertical line,passing through the center 
upon the value for the spray 
v~locity. However, 
of the burst :nd a line passing through 
for bursts on the bottom. the yield should be 
the point. of interest on the surface and 
doubled (as an upper limit estimate) when calcu- 
the center of the burst (degrees), 
~pray velocities. 
. 


, " ... Higher velocities than given by the above 
db 
= .the burst depth (ft). 
~" 
' 
equatIon have been observed for deep bursts,· 
therefore a 50 percent safety factor should be 
added to calculate values of initial dome ve- 


~ 4ocit ., 
, 
h = maximum spray dome height (ft). 
An upper limit of dome height h can be 
lve 
from: 
. 
. 
'lk' 
. 
'. 
r 
__ 
- 
_ 
sd 
- 
- 
_ 
h 
I 
- 
h = V t - 
(j)t 2 
, 
0 


where 


t = time (sec), 


f = retardation factor (ft/sec2 ), 


.ChOPPY sea, together with "Taylor instability" a,ppear to 
account qualitati~ly for this observalion. 


I 


~, 
I 
I 
1 


11 .• h 
"h' 
I 
d' 
.~ uy assummg ,t at graVIty IS t e on y retar mg 


force and, therefore, that f = 16 ft/Sec2 , The 
time that the upper limit height is reached can 
be determined that the flTSt derivative of h with 
respect to time (that is, dome velocity) = zero at 
that time, Therefore, 


- 2 (f)t 


,n. 
'-.' 


e III ility: Spray dome angles are 
estimate 
to be reliable to within ±]0 percent 
for db /W1/3 < ],260 ft/(kt)lf3. Initial vertical 
velocities that are calculated by the equations 
given above might be low by as much as SO per- 
cent, therefore, a safety factor is included in the 
method. The retardation factor f is taken to be 
]6 ft/sec2 to obtain an upper limit for dome 
.height. Actual retardation factors aslarge as 30 
ft/ser? have been observed during nuclear tests. 
.• RelDted MateriDl: See paragraph 2-73. 
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Figure 2-122. _, Velocity of Shock Wave in Free Water 
vs Peak Shock Pressure _ 


Problem 2·33. 
Calculation of Plume, Column, 
and Cauliflower Formation 


700 WI/4 = 1,245 ft. 


Since 240 WI / 4 < db < 700 W1/4 , this is a deep 


UUI:>L. ~ ... ,I<.; ......,.:;ci time is 


. _ 
Figure 2-123 shows reduced plume di- 
mensIOns as a function of reduced time for very 
deep and deep explosions. Figure 2-124 shows 
reduCt'd plume dimensions for shallow bursts 
and cloud di.memions. The criteria for detennin- 
ing the explosion categories as a function of 
depth of burst and weapon yield are given in 
Figure 2·105 and Table 2-8. The accuracy of the 
predictions obtained from 
Figures 2·123 and 
2-124 has not been established fl!11lly because of 
the limited num~er of nuclear tests upon which 
Wllf 


d. 
_. 


. 
Example 


n'elL 'A 10 kt explosion at a depth of 500 


ftin deepwater. 
Find: The plume radius and plume height 
IS seconds after the explosion. 
Solution: 


7 sec/(kt)1l3. 
IS 
=_Ill:: 


(l0)1/3 . 


t 


Wl/3 


Plume height 
Ill:: 800 WI /3 ::::. (800)(10)1/3 


== 1,700 ft, 


From Figure 2-] 23, the reduced plume height is 
about 800 ft/(kt)I/3. and the reduced plume 
radius IS about 700ft/(kt)1 /3. 


. Ans~er: 


Plume radius 
::::. 700 WI /3 
::::. (700)( 10)1/3 


Ill:: 1,500 ft.. 
1II1II Reliability: The curves in Figures 2-123 
and 2-124 are estimated to be reliable to within 
±30 percent. For safety considerations, multiply 
the radius and the height by 1.3. For assurance 
of damage, multiply the values obtained from 
tMe 
fj ures by 0.7. 
Related Material: See "SURFACE EF- 
FE 
OTHER THAN WAYES," in particular 
paragraph 2-74. See also Table 2·8 and Figure 
2-105. 
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Figure 2-123.II1II 
Plume Dimensions'for Deep and Very Deep Unde~water Bursts 
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F"igure 2-124_ Column and Cloud Dimensions for 
Shallow and Very Shallow Bursu 
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Problem 2-34. 
Calculation of Base Surge Growth 
and Cloud Height 


The reduced oiJtside radius of vari- 
ous base surge clouds is shown in Figures 2-126 
and 
~-1;,7 as a function of reduced time after 
burst. Figure ;'-1;'6 shows a curve that is applic- 
able to shallow and very shallow bursts, while 


Figure~-I 27 shows separate curves for deep and 


very deep bursts (see Table 2-8 and Figure 2-105 
for burst depth categories). Equations that relate 


the reduced radii and times to actual radii and 
times for various explosion conditions are shown 
below: 


2-126 


2-125 


2-125 


Applicable 
Figure 


Reduced 
Reduced 
Equation for 
Radius, Rsr 
Time, tf 
Estimated" 
(DimensIonless) 
(seclft!.'- ) 
Dmax (ft) or Amax 


R/Dmax 
t (0 
)1,'2 
Dmax - 710 WI!3 .. 
5 
max 


Rs/Dmax 
t (D 
)1/2 
Dmax 
377 WI !3 (db /wI/3)1/6 
5 
max 


'Rs/Ama\ 
t (A 
. )112 
Amax 
1,50.0 (Wlll l/3 


5 
max 


Burst 


ClassifJca tlon 


Very shallow 


Shallow 


Deep and very 
.deep 


• 
. 
DThe equation D 
ax =710 WI !3 IS valid for very shallow explosions off the bott~m. The equation 
D 
. =377 WI/3 (db/W1fJ\1/6 is villid for very shallow and shallow explosions on a bottom, 
max, 
, 


ZI~ 
2-'2!5- . 


75 WI/3 = (75)(40)1/3 = 256 ft,' 


W = yield of the weapon (kt) 


Rs = radiu~ of the base surge (ft) 


I 
Is = time after explosion (sec) 


~ 
". depth of the explosion (ft) 
-b 


therefore, a 40 kt explosion ata depth of 200 
feet is a very shallow explosion. From the rela- 
tions given above (since the burst was not de- 
scribed as being on the bottom) 


Dmax = 710 WI/3 = (710)(40)1/3 = 2,430 ft, 


Z 
= dbi + 33 =total hydrostatic pressure at 
and 


db: (ft of water) 


= estimated maximum column radius (ft) 
600 
= --:..;.~- = 12.2 sec. 
(2,430)1/2 


I 
f 
1 
,,f 


.. 
_._~max = 'estimated maximum bubble radius (ft). 


_ 
Fig~re 2-125 ~as developed to approxi- 


, 
mate the height of the visible base surge for use 
with radiation exposure calculations. The curve 
in Figure 2·125 was drawn under the assumption 
that the height is not a function .of yield or 
depth of burst. This assumption probably is not 
correct, but there is not sufficient infonnation 


~atethese parameters. 
~Example (U): 
Given: A 40 kt explosion at a depth of 200 
leet. 
Find:, 
a. 
The base surge radius Rs' 
b. 
The height of the visible cloud h 
10 
, 
c 
minutes after burst. 


Solution: Refer to Table 2-8: 


21 WI/3 = (21)(40)1/3 = 72 ft, 


From 
Figure 
2·126, the 
reduced base surge 
radius is about 7. 
Answer: 
a. 
Rs 
=::: Rsr x Dmax = (7)(2,430) 


R 
=::: 17,000 ft. 


i 
' 


b. 
From Figure 2-125, the height of the 
base surge cloud is 


h 
=::: 2,000 ft. 
c 
.... Reliability: 
Base 
surge 
radii 
obtained 


~igures 2·125 and 2-126 are estimated to 
be reliable to within ±20 percent. No estimate 
can be made of the reliability of the height of 
the 
base 
surge 
cloud 
obtained 
from 
Figure 


• 
Related Material: 
See paragraph 2·75. 
See also Table 2-8 and Figure' 2·105. 
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Figure 2-125.11 
Height of Visible Base Surge Cloud 


as a Function of Time ,After Burst., 
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Figure 2-126._ Reduced Base Surge Radius vs Reduced Time 
for Shallow and Very Shallow Bursts • 
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Figure 2~127._" Reduced Base Surge Radius vs Reduced Time 
for Deep and Very Deep Bursls • 
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Problem 2-35. 
Calculation of the Properties of 
Water Surface Waves 


f 
i 


-igUre 2-1.28 shows the peak wave 
~to trough) for any depth of burst in 
shallow water (for waves, shallow water is 0 ..- 
d 91 < 100 W1 / 4 , see paragraph 2-77), and for 
any depth of burst not greater than five times 
the lower cri,tical depth (0'" db /5 <; 170 W1l4 , 
see paragraph 2-77) in deep water (for wavesd91 
~ 256 WI /4 ). The dashed transitional curve in 
Figure 
2-128 
can 
be 
used 
for 
intermediate 
depths. For any depth of burst in shallow water, 
and for bursts as low as the lower critical depth 
(db = 170 WI / 4 ) in deep water, the foregoing 
limits on water depth are also the estimated 


-cavity radii as shown.in Figure 2-128. Fornu- 
clear explosions in deep water at shallow depths 
of burst near the upper critical depth (db 
ll::: 8.6 
Wl l 3 ), the peak wave height may be nearly 
twice as large as predicted in Figure 2-128 (see 
Fi ure 2-107). 
'Exam le 


RelJdblJlty, 


Deep water. 
a. 
The cavity radius giYenby the equation 


, 
',~ , . 


"L 
,C 


Related Material: 
See paragraphs 2·67 
-77. See also Tables 2-9 and 2-10. 


2. Shallow water. 
a. 
The cavity radius given by the equation 
in paragraph 2·77 (and Figure 2-128) provides a 
lower limit for explo'sions in shallow water. An 
upper limit is probably half the maximum bub- 
ble radius given by the equation in paragraph 
2-67. 


b. 
The 
peak· wave height given by the 
equation 
in 
paragraph 
2-77 
(and 
in 
Figure 
2-128) agrees with the measured values within 
about a factor of 2, which is rou 
I 
the un~r­ 


tainty of the data. 


b. 
The peak wave height given in Figure 
2-128 and the equation in paragraph 2·77 agrees 
with the measured H.E. data within a factor of 
two, and within.ercent of the measured 
data of the two nuc ear ex losio 
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Problem 2·36. 
Calculation of Underwater Crater Dimensions 


• 
Figures 2-129 through 2-131 show the 
apparen t 
crater 
radius, 
the 
apparent·· crater 
depth, and the average height of the crater lip, 
respectively, as a function of weapon yield. Fig- 
ures 2-129a, 2·130a, and 2-131 a are for surface 
bursts. Figures 2-129b, 2-130b, and 2·131b are 
for bursts near a clayey sand bottom. Linear 
interpolation between the curves may be used 


_ 
burs 
between the surface and the bottom. 
The curves in Figures 2-129 through 
2·13 
.were developed from 
nuclear cratering 
data; however, H.E. data were used to obtain 
scaling laws and the shape of the curves for bot- 
tom placement. Data from about 600 H.E. and 
.uclear explosions were used in the prepara- 
tion of the curves. The nuclear data are shown in 
Table. 2-11, 


The critical water depth dividing washed 
nwashed craters can be found from Figure 
2-132 for a clayey-sand bottom. An unwashed 
crater results when the water depth is less than 
the minimum height of the irregular cavity lip 
and when the lip does not fail. For other bottom 
materials, the critical water depth can be found 
by multiplying the value of de by the appro- 
priate soil factor for crater lip height in Table 


~Example _ 
~. A 40 kt ~sion on the bottom, 
which consists of oceanic ooze, in 60 feet of 
water. 
Find: The radius, depth, and lip height of 
the apparent crater. 


Solut;on.~ From Figures 2-129,2·130, and 
2-13 I, the crater dimensions in clayey sand are 
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Table 2-12.• SOil' Factors for Cratering in 
Various Materials as COiiared' 


with Clayey Sand 


dae 
~ (0.2)(42) 
~ 8 ft, 


hR, 
~(O.I)(1.7) 
~ 0.2. ft. 


+1500/~50% 
.radius 


Dimensions from these curves are con- 
sidered to have the following reliabilities: 


radius, rae 


~ctors estimated or Interpolated from the 
others based on experimental data.. 


Soil Factors 


Crater 
Crater 
Crater Lip- 
Material. 
Radius 
Depth 
Height· 


Clayey san-d 
1.0 
1.0 
1.0 


Sand 
l.l 
0.8 
0.7 


Loess 
0.9 
1.3 
1.2 


Muck" 
1.5 
0.5 
0.4 


Clayey silt 
].1 
1.1 
1.0 


Coral sand 
0.9 
0.8 
0.7 


Cia: 
0.9 
1.4 
1.3 


Soft rock 
0.8 
0.6' 
0.6 


Hard rock" 
0.6 
0.4 
0.4 . 


Oceanic ooze" 
1.7 
0.2 
0.1 


depth, dae 
~ 
4~ 
ft, 
depth 
-+ 1500/0-50% 


lip height, h£ 
~ 1.7 ft. 
lip height 
+1600/~60% 


Answer: The corresponding dimensions for 
oceanic ooze are (Table 2·12) 


. . rae 
~ (1.7)(750) 
::::: 1,275 ft, 


The accuracy of the soil factors shown in Table 


~12isunknown. 
. 
• 
Related Material: See "UNDERWATER 
CRATERING," page 2-223. 
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2-130. • 
Apparent Crater Depth as a Function of Yield for Various Water Depths 
Over a Clayey Sand Bonom, where dw = Water Depth in Feet. 
Shad Areas 
Denote Possible Transition Regions from a Washed to an Unwashed 
Crater for dw = 20, 40, and 60. 
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